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ABSTRACT 

A  SMALL  COMPUTER  EXPERT  SYSTEM  FOR  LOW-LEVEL  TURBULENCE 
FORECASTS  AT  FORT  IRWIN,  CALIFORNIA 

by  Nelson  L.  Smith 

This  thesis  addresses  the  need  for  comprehensive, 
accurate  Low  Level  Turbulence  (LLT)  forecasts  at  the  US 
Army's  National  Training  Center  (NTC) ,  Ft  Irwin,  California. 
A  LLT  forecast  expert  system  was  developed  for  use  on  a 
small  computer.  The  program  initially  derives  a  single 
macroscale  turbulence  index  from  observations  of  current 
atmospheric  lapse  rate,  pressure  tendency,  wind  speed,  and 
terrain  roughness  for  NTC.  Subsequently,  a  field  of  values 
for  a  Local  Scale  Turbulence  Index  (LTI)  is  computed  as  a 
function  of  macroscale  index,  local  terrain  roughness  at 
1  km  intervals,  and  winds  generated  from  local  observations 
and  a  mass  consistent  wind  model.  LTI  is  modified  for 
terrain  wake  effects  and  Pilot  Reports  (PIREPS)  of 
turbulence  in  the  area.  Threshold  LTI  values  representing 
turbulence  categories  are  graphically  displayed  on  the 
computer  terminal,  superimposed  on  terrain  contours. 

Current  PIREPS  are  also  displayed.  Verification  data 
indicate  the  model  is  a  useful  operational  tool. 
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1.  Introduction 

1.1  Problem  Definition 

The  impact  of  turbulence  on  aircraft  operations  remains 
one  of  the  most  serious  problems  in  aviation  today. 
Turbulence  accounts  for  24%  of  large  commercial  carrier  ac¬ 
cidents,  and  54%  of  all  weather-related  accidents  (McLean, 
1986) .  It  follows  that  there  is  an  ongoing  need  for  accu¬ 
rate  forecasting  of  turbulence  at  all  levels  in  the  atmo¬ 
sphere,  not  only  for  commercial  and  military  aviation 
(Chandler,  1986;  Miller,  1986)  but  also  for  the  launch  and 
recovery  of  space  vehicles  such  as  the  space  shuttle 
(Kolczynski  et  al.,  1986;  Endlich,  1989). 

Low  level  turbulence  (LLT) ,  which  can  be  defined  simply 
as  bumpiness  in  flight  through  the  planetary  boundary  layer 
(PEL) ,  is  more  commonly  encountered  than  any  of  the  other 
turbulence  types,  not  only  because  all  flights  must  pass 
through  the  PEL,  but  because  the  surface  is  a  source  of  tur¬ 
bulent  eddies,  and  the  primary  response  scales  of  aircraft 
are  in  the  size  range  of  those  eddies.  As  with  CAT,  cases 
of  severe  LLT  (i.e.,  resulting  in  a  momentary  loss  of  con¬ 
trol  of  the  aircraft)  are  infrequent,  however  the  proximity 
of  the  ground  allows  less  room  for  recovery,  making  LLT  a 
more  serious  hazard  than  turbulence  at  higher  levels.  Thus 
accurate  prediction  is  a  high  priority. 
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Despite  this  need,  current  LLT  forecasting  methods,  es¬ 
pecially  for  dry  convection  and  mechanical  turbulence,  are 
crude  and  subjective  at  best.  McLean  (1986)  has  pointed  out 
that  although  the  physical  causes  of  LLT  are  generally  well 
understood,  forecast  errors  still  result  from  several 
causes,  including  differences  in  scale  between  observations 
and  forecast  parameters,  initial  values  which  themselves 
must  be  estimated,  limited  data  for  verification,  differing 
levels  of  forecaster  experience,  and  subjective  forecast 
terminology.  These  problems  are  exacerbated  in  the  presence 
of  complex  terrain. 

An  example  of  a  specif' c  LLT  forecast  problem  is  found 
at  the  US  Army  National  Training  Center  (NTC)  at  Ft  Irwin, 
CA.  NTC  is  charged  with  providing  realistic  combat  training 
to  Army  units,  emphasizing  coordinated  air  and  ground  tac¬ 
tics.  An  average  of  3000  sorties  per  year  (half  of  all  tac¬ 
tical  Army  air  operations  in  the  US)  are  flown  at  NTC.  Most 
air  operations  are  conducted  in  helicopters  below  50  feet 
AGL,  and  LLT  prediction  is  a  major  concern. 

Fort  Irwin  is  a  2600  km^  area  in  the  Mojave  Desert  of 
California  (Figure  1) .  It  is  characterized  by  a  variety  of 
terrains  including  high  mountain  peaks,  broad  valleys,  and 
dry  lake  beds  (Figures  2  and  3) .  Because  of  its  location, 
NTC  experiences  a  variety  of  weather  and  associated  LLT:  in¬ 
tense  dry  convective  activity,  thunderstorms,  strong  winds 
associated  with  frontal  passages,  and  mountain  waves.  LLT 


Figure  1  Ft  Irwin  Location 

Circles  represent  surface  reporting  stations.  Indicating 

call  sign  and  number. 


Altitude 
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Ft  Irvnn  Terrain  (200  m  contours) 


Figure  2  Three-Dimensional 
View  of  Ft  Irwin  Terrain 

Vertical  scale  exaggerated  10  times. 

west . 


View  is  from  south- 


Figure  3  Ft  Irwin  Topography 


Contours  are  at  200  meter  Intervals,  axes  are  In  Unlve 
Transverse  Mercator  CUTM5  kilometer  coordinates 
Meteorological  sensors  are  located  on  Towers  1-5  . 
Indicated.  BYS  Is  Bicycle  Lake  Army  Air  Field. 


forecasts  are  complicated  by  lack  of  data  and  the  rugged 
terrain  prevalent  at  NTC.  (Air  Force,  1987) 

Point  warnings  for  NTC  are  issued  by  the  US  Air  Force 
Global  Weather  Central  (GWC) ,  with  supplemental  forecasts 
issued  by  George  AFB  (VCV) ;  operational  forecasts  are  pro¬ 
vided  by  forecasters  attached  to  (transient)  trainee  units. 
The  lack  of  forecaster  experience  compounds  the  LLT  forecast 
problem  at  NTC.  Transient  forecasters  have  approximately 
three  weeks  to  become  familiar  with  the  terrain  and  weather 
patterns  common  at  NTC.  Their  natural  tendency  is  to  fore¬ 
cast  conservatively  in  unfamiliar  situations;  this,  coupled 
with  an  inability  to  resolve  the  scales  involved  and  a  lack 
of  local  observations,  results  in  the  user's  perception  that 
the  forecasts  are  inaccurate  (i.e.  over-forecasts)  (Lester 
and  Burton,  1988) .  For  example,  an  average  of  20%  of  all 
missions  during  spring  and  fall  are  cancelled  due  to  fore¬ 
cast  LLT  and/or  predicted  high  surface  winds  (Marrs,  1988) . 
Although  some  of  these  cancellations  are  clearly  warranted, 
in  many  cases  only  a  portion  of  NTC  experiences  LLT  of  crit¬ 
ical  intensities;  i.e.,  flying  activities  could  be  conducted 
safely  in  other  areas.  Current  LLT  forecasts  do  not  contain 
such  detail  and  the  post  is  either  entirely  open  or  closed 
to  flying. 

Recently,  several  experiments  were  carried  out  in  an 
attempt  to  overcome  the  LLT  forecast  problem  at  NTC:  (i)  a 
network  of  automated  surface  weather  stations  was  installed 
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(Marrs,  1988) ;  (ii)  a  small  computer  scheme  for  real-time 
objective  analysis  and  graphic  display  of  terrain  and  winds 
across  the  post  was  developed  (Henmi  et  al.,  1987);  and 
(iii)  an  automated  (small  computer  based)  system  for  general 
(point)  forecasts  of  LLT  at  Ft  Irwin  was  developed  (Lee, 
1988).  In  addition,  Lester  and  Burton  (1988):  (iv)  mapped 
turbulence  prone  areas  at  NTC  (Appendix  I) ;  (v)  initiated  a 
turbulence  data  collection  scheme;  and  (vi)  expanded  on  ear¬ 
lier  work  by  Burton  (1964)  and  Ludwig  and  Endlich  (1988)  to 
demonstrate  the  feasibility  of  the  real-time  computation  and 
display  of  a  map  of  an  objective  index  of  LLT  intensity  for 
the  local  area. 

These  initial  experiments  have  met  with  varying  degrees 
of  success.  For  example: 

The  increased  number  of  observation  stations,  and  the 
real-time  display  of  the  data  on  a  map  of  the  post  was  found 
exceptionally  useful  by  the  forecasters  and  pilots  alike. 
However,  better  coverage  of  critical  areas  is  needed. 
Although  the  interpolation  of  data  from  the  improved  network 
to  all  locations  across  the  post  via  the  Henmi  (1987)  mass 
consistent  wind  model  provided  the  needed  coverage,  the  in¬ 
fluence  of  stability  on  the  distribution  of  winds  was  ig¬ 
nored  in  the  model. 

The  LLT  forecast  aid  developed  by  Lee  (1988)  provides  a 
useful  training  tool  for  forecasters  unfamiliar  with  the 
area,  as  well  as  a  semi-objective  forecast  "checklist”  which 
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helps  to  standardize  forecast  methods  and  procedures.  The 
main  disadvantage  of  this  method  is  that  it  yields  only 
gross  point  forecasts  which  do  not  address  the  problem  of 
turbulence  gradients  across  the  post. 

The  turbulence  climatology  developed  by  Lester  and 
Burton  (1988)  is  a  useful  familiarization  tool  for  new  fore¬ 
casters.  Also,  their  turbulence  data  collection  program  has 
resulted  in  a  useful  data  set  for  winter  and  spring,  1988, 
at  NTC.  Their  turbulence  index  scheme,  however,  has  not 
been  tested  against  extensive  data. 

The  results  of  these  preliminary  experiments  present  a 
starting  point  for  the  significant  improvement  of  LLT  fore¬ 
casts.  The  ability  to  manipulate  large  data  bases  and  to 
execute  relatively  sophisticated  physical  models  locally  on 
small  computers  allows  the  forecaster  to  make  more  detailed 
analyses  than  are  practical  at  a  central  facility. 
Furthermore,  the  incorporation  of  artificial-intelligence 
(AI)  technology  into  forecast  schemes  offers  the  potential 
for  a  marked  improvement  in  forecast  accuracy  and  consis¬ 
tency  (Racer  and  Gaffney,  1986) . 

1.2  Objective 

The  previous  section  has  documented  the  need  for  re¬ 
search  toward  better  LLT  forecasts.  Furthermore,  the  focus 
of  recent  LLT  research  at  Ft  Irwin  has  demonstrated  improve- 
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merits  which  may  be  expected  with  the  application  of  avail¬ 
able  small  computer  technology.  The  objective  of  the 
current  research  is  to  integrate  and  expand  the  NTC  work  to 
develop  a  small  computer  expert  system  for  LLT  forecasts. 

1 . 3  Scope 

The  proposed  expert  system  will  be  based  on  the  objec¬ 
tive  low  level  Local-Scale  Turbulence  Index  (LTI)  proposed 
by  Lester  and  Burton  (1988).  A  brief  review  of  current  LLT 
forecast  methods,  the  prototype  LTI  system,  and  AI  methods 
will  be  presented  in  Chapter  2.  System  development  is 
described  in  Chapter  3,  with  a  description  of  available  test 
data  in  Chapter  4.  Program  execution  and  results  are 
presented  in  Chapter  5  followed  by  conclusions  and 
recommendations  (Chapter  6) . 
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2 .  Background 

2 . 1  Low  Level  Turbulence 

’’Turbulence”  is  commonly  defined  in  terms  of  fluid  be¬ 
havior.  For  example,  Huschke  (1959)  defines  turbulence  as 
”a  state  of  fluid  flow  in  which  the  instantaneous  velocities 
exhibit  irregular  and  apparently  random  fluctuations  so  that 
in  practice  only  statistical  properties  can  be  recognized 
and  subjected  to  analysis.”  In  contrast,  aviation  turbu¬ 
lence  is  defined  with  respect  to  the  aircraft  which  encoun¬ 
ters  it,  e.g.,  as  "bumpiness  in  flight.”  It  follows  that 
the  intensity  of  aviation  turbulence  is  usually  expressed  in 
terms  of  its  effects  on  aircraft  (Table  1) .  In  fact, 
semi-quantitative  £ilot  Reports  (PiREPs)  of  those  effects 
are  the  most  common  form  of  aviation  turbulence  data 
available  for  both  forecasts  and  research. 

Low  Level  Turbulence  (LLT)  or  "bumpiness  in  flight 
through  the  planetary  boundary  layer"  is  caused  by  a  number 
of  phenomena,  including  dry  convection  (thermals) ,  moist 
convection  (thunderstorms,  downbursts,  etc.),  mechanical 
mixing,  mountain  waves,  low-level  wind  shear,  and  fronts. 
These  mechanisms  are  not  mutually  exclusive,  but  their  sepa¬ 
ration  is  useful  in  the  discussion  of  turbulent  processes. 
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Table  1  Turbulence  Categories  (AWS,  1988) 


Light  turbulence  may  cause  slight  erratic  changes  in  atti¬ 
tude  and/or  altitude.  It  produces  a  variation  in  airspeed 
from  5  to  15  knots.  Seat  belts  may  be  required  and  occu¬ 
pants  may  feel  a  slight  strain  against  restraints.  Loose 
objects  in  the  aircraft  may  be  displaced  slightly.  Food 
service  may  be  conducted  and  little  or  no  difficulty  is  en¬ 
countered  while  walking. 

Moderate  turbulence  causes  changes  in  altitude  and/or  atti¬ 
tude  but  aircraft  remains  in  positive  control.  Airspeed  is 
affected,  varying  from  15  to  25  knots.  Occupants  feel  def¬ 
inite  strain  against  restraints.  Unsecured  objects  are  dis¬ 
lodged.  Food  service  and  walking  are  difficult. 

Severe  turbulence  causes  large,  abrupt  changes  in  altitude 
and/or  attitude.  Airspeed  is  affected  in  excess  of  25 
knots.  Aircraft  may  be  momentarily  out  of  control. 

Occupants  are  forced  violently  against  seat  belts  and  the 
seat.  Loose  objects  are  tossed  about  the  aircraft.  Food 
service  and  walking  are  impossible. 

Extreme  turbulence  is  very  rare.  The  aircraft  is  violently 
tossed  about  and  is  practically  impossible  to  control. 

Rapid  fluctuation  of  airspeed  occurs  in  excess  of  25  knots. 
It  may  cause  structural  damage. 
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2.1.1  LLT  Forecast  Methods 

Current  LLT  forecast  procedures  rely  on  the  association 
of  one  or  more  large  scale  meteorological  indicators  with 
the  occurrence  of  the  mechanisms  listed  above.  Examples  are 
large  scale  flow  patterns  (FAA,  1987;  AWS,  1988;  Lee,  1988), 
calculated  parameters  or  indices  (Burnett,  1970;  Lester  and 
Burton,  1988) ,  or  rules  of  thumb  often  related  to  a  thresh¬ 
old  value  of  some  meteorological  variable  (e.g.  Lee  et  al., 
1979)  . 

Recognizing  the  subjectivity  inherent  in  the  applica¬ 
tion  of  these  forecast  procedures,  Lester  and  Burton  (1988) 
have  presented  a  simplified  flow  diagram  (Figure  4)  to 
clarify  the  mental  processes  meteorologists  use  to  formulate 
LLT  forecasts.  This  procedure,  although  idealized,  is 
comprehensive  and  will  serve  as  a  framework  for  further 
discussion. 

As  shown  in  Figure  4,  "Pattern  Recognition”  refers  to 
the  identification  of  synoptic  scale  (and,  where  possible, 
mesoscale)  flow  patterns  associated  with  LLT.  The  synoptic 
scale  circulation  patterns  associated  with  LLT-generating 
phenomena  have  been  well  described.  Synoptic  patterns  fa¬ 
vorable  for  dry  convection  have  been  examined  in  the  context 
of  forecasts  for  soaring  (Wallington,  1966;  Lindsay  and 
Lacy,  1976;  Bradbury  and  Kuettner,  1976).  Conditions  con¬ 
ducive  to  the  development  of  moist  convection  have  long  been 


Figure  4  Idealized  LLT  Forecast  Procedure 
After  Lester  and  Burton  C1988) 
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of  interest  to  operational  forecasters,  especially  with 
respect  to  severe  weather,  and  have  been  described 
repeatedly  in  the  literature  (e.g..  Miller,  1972;  Ray,  1986; 
AWS  1988) . 

Waters  (1970) ,  among  others,  has  reviewed  synoptic 
patterns  frequently  associated  with  strong,  gusty  winds  of 
non-convective  origin.  Mountain  wave  patterns  have  been 
studied  extensively  with  reviews  by  Alaka  (1960) ,  Nicholls 
(1973) ,  and  others.  Frontal  patterns  are  also  described  in 
numerous  sources  such  as  Pettersen  (1956) ,  Palmen  and  Newton 
(1969),  Byers  (1974),  Chandler  (1986),  and  AWS  (1988). 
Additionally,  local  Terminal  Forecast  Reference  Notebooks 
(TFRN)  typically  include  local  patterns  favorable  to  the 
development  of  one  or  more  of  the  conditions  listed  above. 

With  the  exception  of  radar  and  satellite  observation 
of  moist  convection,  mesoscale  circulation  patterns  are  not 
as  well  measured  as  synoptic  patterns.  Lester  and  Burton 
(1988)  note  that,  to  fill  this  gap,  an  individual  forecaster 
will  develop  a  "Conceptual  Model"  based  on  accepted  theory, 
on  experimentally  developed  structural  models,  or  simply  on 
experience  (Figure  4) . 

Mesoscale  models  of  dry  convection  are  described  in  the 
gliding  literature  cited  earlier.  Moist  convective  models 
have  been  well  described  by  Palmen  and  Newton  (1969) ,  Byers 
(1974),  Atkinson  (1981),  Fujita  (1985),  Ray  (1986)  and 


others. 
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Structural  models  of  mountain  wave  systems  have  been 
applied  extensively  to  infer  areas  of  turbulence  (Alaka, 
1960;  Nicholls,  1973;  Lester  and  Fingerhut,  1974),  and 
severe  downslope  windstorms  (Lilly,  1978;  Durran,  1986; 
Giusti,  1987)  for  many  years. 

Useful  guidance  for  anticipating  the  location  and 
intensity  of  mechanically  generated  turbulence  and 
turbulence  in  the  vicinity  of  fronts  has  received  less 
attention  (Theon,  1986;  Chandler,  1986;  AWS,  1988) . 

Once  the  likelihood  of  LLT  is  recognized,  some  sort  of 
"Parameter  Evaluation”  (Figure  4)  is  used  to  determine  the 
exact  areas,  times,  and  intensities  of  the  expected 
turbulence.  Lester  and  Burton  (1988)  group  these  parameters 
into  three  categories:  (i)  the  basic  meteorological 
variables,  (ii)  their  spatial  and  temporal  derivatives,  and 
(iii)  physical  and/or  empirical  indices. 

The  forecast  parameters  associated  with  moist 
convection  are  extensive  and  well  known  (see,  e.g.,  Ray, 
1986;  AWS,  1988).  They  include  numerous  stability  indices 
and  combinations  of  basic  variables  such  as  temperature, 
relative  humidity,  and  wind  shear.  In  addition, 
observational  tools  such  as  radar  and  satellites  are  useful 
in  determining  the  location  and  intensities  of  thunderstorm 
elements. 

Table  2  lists  some  common  forecast  tools  for  LLT  not 
associated  with  thunderstorms.  The  simpler  parameters 
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Table  2  Some  Common  LLT  Forecast  Tools 
(After  Lester  and  Burton,  1988) 

Dry  Convection  fTherroall  LLT 


Surface  Temperature 
Temperature  Lapse  Rate 
Potential  Temperature  Lapse  Rate 
Thermal  Index  (Higgins,  1963) 


rate 

Mechanical  LLT 

Surface  Wind  Speed  and  Gusts 
Gradient  Level  Wind  Speed 
Terrain  Roughness 
Mountain  Wave  LLT 


T 

■^0 

6T/OZ 

66/6z 

T,  -  (Tq  +  rz) 

T^  =  Temperature 
at  850  mb 

r  =  Dry  adiabatic 
lapse 


V, 


V 


Mountain  Top  Wind  Speed 

Cross  Mountain  Sea  Level  Pressure 

Gradient  6P 

LLT  Incites 


Panofsky  Index  PI  =  V  (l  -  Ri/Ri^.^) 

(US  Navy,  1975) 

V  =  mean  wind 
speed 

Ri  =  gradient 

Richardson  # 
for 

boundary  layer 

AFGWC  Mechanical  Turbulence  Index  I  =  aUR  +  b 

(Burnett,  1970) 

a,  b  =  constants 
U  =  mountain  top 
wind  speed 
R  =  upwind  terrain 
roughness 
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(e.g.,  surface  temperatures)  are  often  used  because  the 
vertical  structure  of  the  atmosphere  in  the  area  of  interest 
is  rarely  known  in  great  detail.  The  combination  of 
parameters  is  frequently  accomplished  via  indices  (such  as 
those  shown  at  the  bottom  of  Table  2) ,  via  "look-up”  tables 
(e.g.,  Table  3),  or  with  nomograms.  Parameters  for 
estimation  of  LLT  in  the  vicinity  of  large-scale  fronts 
usually  depend  on  some  knowledge  of  the  intensity  and  speed 
of  the  front  (e.g.,  Figure  5).  Richwien  (1979)  has  shown  an 
association  between  significant  LLT  and  frontal  systems  with 
horizontal  temperature  differences  of  at  least  10  degrees  F 
and  frontal  speeds  of  at  least  30  knots.  Fronts  moving 
across  rough  terrain  further  enhance  production  (Lester 
and  Burton,  1988) . 

Threshold  values  for  critical  parameters  are  not 
universal;  they  are  usually  poorly  defined  functions  of 
location,  season,  time  of  day,  and  aircraft  type.  Lee  et 
al.  (1979),  the  FAA  (1987)  and  the  AWS  (1988)  give  typical 
values  for  these.  The  FAA  (1987)  also  gives  a  useful 
summary  of  general  rules  of  thumb  derived  from  these 
parameters  (Table  4) .  In  addition,  Lester  and  Burton  (1988) 
note  two  other  well-known,  but  frequently  unstated,  rules  of 
thumb;  (i)  the  intensity  of  turbulence  always  increases  in 
proportion  to  wind  speed  and  roughness;  and,  (ii)  turbulence 
elements  (hence  surface  gusts  and  LLT)  have  dimensions 
proportional  to  the  size  of  the  roughness  elements. 
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Table  3  Low  Level  Turbulence  Table 
AWS  C1988D 


Stability 


Surface  Wind  Speed 


1-12  13-24  25-50  >  50 

kts  ktg  kta  kta 


Smooth  Terrain  Very  aeable  0  I.  L  H 

Relatively  (table  0  L  H  M 

Relatively  unstable  I.  H  S  S 

Very  unstable  H  M  S  S 


Roush  Terrain  Very  (table  0  L  M  S 

Relatively  (table  0  H  S  S 

Relatively  unstable  L  K  S  S 

Very  unstable  H  S  S  S 


Letters  Indicate  turbulence  intensity  as  follows: 


0;  None  L:  Light  M:  Moderate  S:  Severe 

Rough  terrain  1(  defined  as  a  region  where  topogrephle  features  extend  sore  than  1,000  feet  above  the 
surroundings. 

The  wind  velocity  used  Is  the  surface  wind  Including  gusts. 

The  stability  categories  apply  to  the  layer  from  the  surface  to  850  mb  and  are  defined  as  follows: 

Very  stable:  A  lapse  rate  equal  to  or  less  than  the  moist  adiabatic  race. 

Relatively  stable:  A  lapse  rate  between  Che  moist  adiabatic  and  the  mean  lapse  race.  The 
mean  lapse  rate  la  defined  as  that  midway  between  the  moist  and  dry  adiabatic  races. 

Relatively  unstable:  A  lapse  race  between  the  mean  and  dry  adiabatic  races. 

Very  unstable:  Lapse  race  equal  to  or  greater  then  the  dry  adiabatic  race. 


Temperature  difference  acroa 
front  at  surface  (degrees  F) 

(AT-  ®F/50NM) 
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TURBULENCE  ASSOCIATED  WITH  FRONTAL  LOW  LEVEL  WIND  SHEAR 


Figure  5  Turbulence  Associated  With  Frontal  Low  Level  Wind 

Shear 


AWS  C1987D 


Table  4  Locations  of  Probable  Turbulence  by  Intensities 
Versus  Weather  and  Terrain  Features  (FAA,  1987) 
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Light  Turbulence 

1.  In  hilly  and  mountainous  areas  even  with  light  winds. 

2.  In  and  near  small  cumulus  clouds. 

3.  In  clear-air  convective  currents  over  heated  surfaces. 

4.  In  the  lower  5,000  ft  of  the  atmosphere: ^ 

a.  When  winds  are  near  15  kts. 

b.  Where  the  air  is  colder  than  the  underlying 
surfaces. 

Moderate  Turbulence 

1.  In  mountainous  areas  with  a  wind  component  of  25  to  50 
kts  perpendicular  to  and  near  the  level  of  the  ridge: 

a.  At  all  levels  from  the  surface  to  5000  ft  above 
the  tropopause  with  preference  for  altitudes: 

(1)  Within  5000  ft  of  the  ridge  level. 

(2)  At  the  base  of  relatively  stable  layers  below 
the  base  of  the  tropopause. 

(3)  Within  the  tropopause  layer. 

b.  Extending  outward  on  the  lee  of  the  ridge  for  150 
to  300  miles. 

2.  In  and  near  thunderstorms  in  the  dissipating  stage. 

3.  In  and  near  other  towering  cumuliform  clouds. 

4.  In  the  lower  5,000  ft  of  the  troposphere: 

a.  When  surface  winds  exceed  25  kts. 

b.  Where  heating  of  the  underlying  surface  is 
unusually  strong. 

c.  Where  there  is  and  invasion  of  very  cold  air. 

5.  In  fronts  aloft. 

6 .  Where ; 

a.  Vertical  wind  shears  exceed  6  kts/1000  ft,  and/or 

b.  Horizontal  wind  shears  exceed  18  kts/150  miles. 
Severe  Turbulence 

1.  In  mountainous  areas  with  a  wind  component  exceeding  50 
kts  perpendicular  to  and  near  the  level  of  the  ridge: 

a.  In  5,000-ft  layers: 

(1)  At  and  below  the  ridge  level  in  rotor  clouds 
or  rotor  action. 

b.  Extending  outward  on  the  lee  edge  of  the  ridge  for 
50  to  150  miles. 

Extreme  Turbulence 

1.  In  mountain  wave  situations,  in  and  below  the  level  of 
well-developed  rotor  clouds.  Sometimes  it  extends  to 
the  ground. 

2.  In  severe  thunderstorms  (most  frequently  in  organized 
sc[uall  lines)  indicated  by: 

a.  Large  hailstones  (3/4  inch  or  more  in  diameter) . 

b.  Strong  radar  echoes,  or 

c.  Almost  continuous  lightning. 
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Refinement  of  these  rules  of  thumb  for  specific  areas 
occurs  during  the  "Local  Tuning"  process  (Figure  4) ,  which 
is  based  on  verifications  from  occasional  pilot  reports 
(PIREPS) ,  knowledge  of  local  conditions,  and  experience. 

Most  rule  modifications  result  simply  in  an  increase  or 
decrease  of  a  critical  threshold  value. 

The  "Metwatch"  function  illustrated  in  Figure  4  entails 
monitoring  conditions  after  a  forecast  has  been  made.  Along 
with  "Local  Tuning,"  "Metwatch"  is  one  of  the  most 
labor-intensive  parts  of  the  overall  forecast  process 
(Richwien,  1979) ;  and,  as  pointed  out  by  Lester  and  Burton 
(1988) ,  the  manpower  required  for  this  function  is  often 
unavailable,  especially  in  critical  forecast  situations. 

The  last  step  in  the  forecast  process  is  comprehensive 
"Verification"  (Figure  4) .  This,  of  course,  is  required  to 
monitor  skill  and  improve  the  quality  of  the  forecast 
product  (McGinley,  1986) .  Due  to  the  nature  of  low  level 
turbulence  data,  this  important  step  in  the  forecast 
procedure  is  often  neglected. 

2.1.2  LLT  Forecast  Problems 

The  primary  problems  associated  with  LLT  forecasting 
are  (i)  scale,  (ii)  forecast  subjectivity,  and  (iii) 
verification. 
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The  critical  time  and  space  scales  of  turbulence  are 
several  orders  of  magnitude  smaller  than  available 
measurements.  For  example,  Roeder  and  Gall  (1987)  have 
estimated  the  width  of  a  typical  cold  front  as  less  than  20 
kilometers.  Jones  et  al.  (1970),  Murrow  (1986),  and  others 
have  shown  the  critical  horizontal  scale  length  of 
atmospheric  turbulence  to  be  on  the  order  of  one  kilometer 
or  less,  and  time  scales  on  the  order  of  minutes.  Terrain 
scales  important  to  turbulence  generation  may  also  be  as 
small  as  ten  meters  (Theon,  1986) .  As  a  result  of 
observational  inadequacies,  important  turbulence  variables 
such  as  atmospheric  lapse  rate  (stability) ,  wind  speeds, 
pressure  and  temperature  gradients  must  be  interpolated  or 
estimated  in  some  way  before  their  inclusion  in  any  LLT 
forecast  scheme.  Clearly,  this  process  decreases  the 
accuracy  of  the  forecast  product. 

Subjectivity  affects  forecasters  and  pilots  alike.  The 
"Pattern  Recognition”  and  Local  Tuning"  tasks  listed  earlier 
are  highly  subjective,  and  can  vary  significantly  with 
forecaster  experience,  training,  and  familiarity  with  the 
local  area.  AWS  (1988)  cites  the  natural  trend  to 
over-forecast  certain  categories  and  under-forecast  others. 
The  turbulence  categories  themselves  (e.g.  "Moderate,” 
"Extreme,”  Table  1)  introduce  further  subjectivity  into  the 
forecast.  Brown  and  Murphy  (1982)  have  shown  the 
superiority  of  a  numerical  index  to  verbal  categories  which 
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are  likely  to  mean  different  things  to  different  people. 
Participants  at  a  recent  turbulence  workshop  conducted  by 
the  National  Aeronautics  and  Space  Administration  (NASA) 
agreed,  listing  the  standardization  and  quantification  of 
turbulence  intensity  terms  as  ’'high  priority"  (NASA,  1986)  . 

"Verification"  is  another  problem  area  in  the 
development  of  an  acceptable  LLT  forecast  technique. 

Accurate  turbulence  data  bases  are  difficult  to  collect, 
given  the  expense  of  instrumented  observations  and  the 
subjectivity  of  PIREPs.  To  date,  the  verification  of  LLT 
forecasts  has  depended  almost  solely  on  the  latter.  Keller 
(1986)  has  emphasized  the  subjective  nature  of  PIREPs, 
pointing  out  that  pilots  aren't  "blind"  sensors,  they  react 
to  encountered  turbulence,  and  avoid  areas  where  it  is 
forecast.  Furthermore,  during  a  turbulence  encounter,  pilot 
corrective  action  can  actually  increase  aircraft  buffeting 
due  to  the  phase  relationship  between  the  turbulence  and 
corrective  action.  These  problems  are  exacerbated  by 
differences  in  pilot  experience,  and  aircraft  type  and 
speed.  Despite  this  subjectivity,  PIREPs  remain  the  most 
effective  means  of  collecting  large  amounts  of  data  over  a 
given  area. 
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2.2  Artificial  Intelligence  (AI)  and  Expert 
Systems  (ES) 

Recent  advances  in  small  computer  technology  provide 
new  opportunities  for  improved  local  forecasts.  Desk  top 
computers  now  have  the  capacity  to  handle  large  data  bases 
and  to  execute  relatively  sophisticated  physical  models  in 
near  real  time.  AI  technology  promises  to  increase  forecast 
consistency  and  reliability  by  bringing  all  forecasters  up 
to  the  same  level  of  experience.  Smith  (1988)  has  recently 
reviewed  AI  applications  to  forecast  problems,  and  the 
following  is  based  partially  on  that  review. 

Artificial  Intelligence  (AI)  is  a  generic  term 
referring  to  the  use  of  a  computer  to  imitate  human  behavior 
which  is  generally  thought  to  require  intelligence.  Expert 
Systems  (ES)  and  Knowledge  Based  Systems  (KBS)  are  less 
stringent  terms  dealing  with  the  use  of  computers  to  emulate 
human  thought  processes  under  stricter  guidelines  using 
empirical  relationships  based  on  experience  and  knowledge  of 
the  programmer.  Racer  and  Gaffney  (1984)  introduce  the  term 
Interpretive  Processing  (IP)  as  an  application  of  ES/KBS  in 
meteorological  applications  and  quote  the  following 
definitions; 

Artificial  Intelligence  is  a  subfield  of  computer 
science  concerned  with  the  concepts  and  methods  of  symbolic 
inference  by  a  computer  and  the  symbolic  representation  of 
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the  knowledge  to  be  used  in  making  inferences.  A  computer 
can  be  made  to  behave  in  ways  that  humans  recognize  as 
"intelligent”  behavior  in  each  other.  (Feigenbaum  and 
McCorduck,  1983) 

According  to  Duda  and  Shortliffe  (1983) ,  Artificial 
Intelligence  is  the  development  of  computer  programs  that 
can  solve  problems  normally  thought  to  require  human 
intelligence. 

Nau  (1983)  defines  Expert  Systems  as  problem-solving 
computer  programs  that  can  reach  a  level  of  performance 
comparable  to  that  of  a  human  expert  in  some  specialized 
problem  domain. 

By  contrast,  Duda  and  Shortliffe  (1983)  define  a 
Knowledge-Based  system  as  an  AI  program  whose  performance 
depends  more  on  the  explicit  presence  of  a  large  body  of 
knowledge  than  on  the  possession  of  ingenious  computational 
procedures;  by  expert  system  they  mean  a  knowledge-based 
system  whose  performance  is  intended  to  rival  that  of  human 
experts . 

Interpretive  Processing  (IP) ,  as  used  here,  is  a 
computer  interactive  procedure  that  enhances  the  abilities 
of  the  weather  forecaster  to  decide  on  a  forecast.  The 
procedure  makes  it  easier  to  draw  conclusions  from  the 
meteorological  analysis  of  observational  data,  forecasting 
techniques,  and  past  forecaster  experience  available  when 
deciding  on  a  forecast. 
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The  possible  applications  of  AI  to  meteorology  cover  a 
spectrum,  ranging  from  decision  trees,  such  as  developed  by 
Brown  (1986)  and  Colquhoun  (1987)  to  forecast  programs 
capable  of  learning  (Gaffney  and  Racer,  1983)  and  beyond. 

The  National  Weather  Service  is  increasing  its  automation  of 
field  operations  as  part  of  its  modernization  efforts,  with 
one  of  its  areas  of  concentration  being  the  field  of 
Interpretive  Processing.  Since  the  forecast  problem 
involves  reduction  of  available  data,  identification  of 
significant  data  and  guidance  (numerical  and  manual) ,  and 
the  application  of  both  explicit  and  implicit  relationships, 
rules  of  thumb,  etc.  to  create  a  forecast  product,  a 
competent  IP  system  would  be  of  great  benefit.  Racer  and 
Gaffney  (1984)  give  an  example  of  a  prototype  IP  system 
tailored  to  NWS  needs.  They  further  envision  a  three-fold 
benefit  from  the  application  of  ES/KBS  to  weather 
forecasting:  (i)  to  provide  improved  data  analysis  and 

decision-making  support  due  to  enhanced  consistency  and 
thoroughness;  (ii)  to  support  training  of  new  forecasters; 
and,  (iii)  to  support  skill  maintenance  for  experienced 
forecasters,  especially  with  regard  to  their  actions  in 
infrequently-occurring,  unfamiliar  situations.  Successful 
incorporation  of  these  objectives  into  a  comprehensive  LLT 
forecast  scheme  would  improve  forecasts  significantly. 

AI  technology  is  being  used  in  varying  degrees  as  a 
forecast  tool.  Brown  (1986)  has  developed  a  simple  decision 
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tree  approach  to  forecasting  downslope  wind  storms  in 
Colorado.  His  is  a  program  using  "if-then”  structures  to 
consolidate  significant  data  (both  analysis  and  numerical 
guidance)  into  a  valid  indicator  of  the  probability  of 
strong  downslope  winds.  Colquhoun  (1987)  has  used  a  similar 
approach  in  the  forecasts  of  thunderstorms  and  tornadoes. 
Gaffney  and  Racer  (1983)  have  developed  a  prototype  system 
for  severe  storm  advisories  which  is  capable  of  "learning" 
behavior.  This  system  is  based  upon  formalized  rules 
developed  by  Crisp  (1959)  and  Miller  (1972)  of  the  Air  Force 
Global  Weather  Central.  Racer  and  Gaffney  (1986)  quote  a 
personal  comm’  -at ion  with  J.T.  Schaefer  of  the  National 

Severe  Stoir -  Forecast  Center  detailing  a  KBS  which  includes 
"a  severe  weather  checklist  of  10  parameters  which  are 
evaluated  as  a  group  using  'if-then'  rules  to  determine  the 
'possibility'  of  a  storm."  Racer  and  Gaffney  (1986)  also 
detail  a  diagnosis  procedure  for  evaluating  numerical 
guidance  materials  developed  by  Simpson  (1971)  at  the  NWS 
National  Hurricane  Center.  It  uses  a  decision  ladder  for 
systematic  analysis  of  the  performance  of  numerical  models 
with  the  goal  of  improving  them. 

There  is  an  apparent  gap  in  the  spectrum  of  technical 
applications  of  AI  to  weather  forecasting.  Gaffney  and 
Racer's  "learning"  program  is  at  the  high  end,  but  it  is 
only  a  prototype.  The  checklist/decision  tree  approach 
(e.g.,  Lee,  1988),  at  the  low  end  of  the  spectrum,  is  the 


only  application  of  AI  currently  in  use.  While  this  is  an 
improvement  over  manual  methods,  much  greater  benefits  could 
be  realized  by  the  use  of  "smarter''  systems. 

2.3  The  Local  Scale  Turbulence  Index  (LTI)  —  A 
Prototype  Expert  System  for  LLT  Forecasts 

The  Burton  Turbulence  Index  is  a  non-dimensional 
numerical  index  used  to  describe  expected  turbulence 
intensities  over  a  large  area.  BTI  uses  three  of  the 
forecast  tools  shown  in  Table  2  to  characterize  LLT:  (i) 
gradient  wind  speed  and  (ii)  surface  roughness  as  indicators 
of  mechanical  LLT,  and  (iii)  temperature  lapse  rate  to 
indicate  dry  convective  or  thermal  LLT.  Moist  convective 
(e.g.,  frontal)  LLT  is  addressed  by  a  fourth  parameter,  the 
surface  pressure  tendency.  BTI  has  been  used  operationally 
as  a  macroscale  LLT  forecast  tool  in  the  past;  Lester  and 
Burton  have  recently  shown  that,  with  modifications,  BTI  may 
be  adapted  to  finer  scales. 

BTI  is  defined  as 

BTI  =  R  +  V  +  S  +  T  (1) 

where  R  is  roughness  (difference  between  the  highest  and 
lowest  elevations)  in  hundreds  of  feet,  V  is  the  wind  speed 
in  knots  at  2000'  AGL,  S  is  ten  times  the  lapse  rate  in 
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•C/1000  ft  in  the  lowest  100  mb,  and  T  is  the  absolute  value 
of  the  3 -hour  pressure  tendency. 

The  BTI  was  used  extensively  from  the  mid  nineteen- 
sixties  through  the  early  seventies  as  a  large  scale  LLT 
indicator  at  the  US  Air  Force  Global  Weather  Central  (GWC) 
(Burton,  1964;  Burnett,  1970).  Table  5  lists  critical  BTI 
values  for  category  I  aircraft  (i.e.,  those  aircraft  most 
susceptible  to  turbulence) . 

Jones  (1970)  has  investigated  the  use  of  wind  speed, 
lapse  rate,  roughness,  BTI,  Richardson  number,  and  Showalter 
index  in  the  prediction  of  LLT  using  data  collected  from  the 
US  Air  Force  LO-LOCAT  project  (Loving,  1969) .  The  results 
clearly  demonstrated  the  usefulness  of  BTI  in  comparison  to 
the  other  variables  and  indices. 

Lester  and  Burton  (1988)  developed  a  prototype 
local-scale  turbulence  index  (LTI)  on  a  mesoscale  (B,t)  grid 
for  NTC.  The  LTI  is  a  scaled  version  of  the  BTI,  and  is 
calculated  as  follows: 

LTI  =  BTI  X  (R.  +  Vp/(R  +  V)^  (2) 

Rj  =  the  roughness  (as  defined  for  BTI)  for  a  two 
kilometer  square  centered  on  grid  point  (j). 

Vj  =  the  10  m  wind  speed  at  grid  point  (j) , 
determined  from  local  observations  and  a 
diagnostic  wind  model. 
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Table  5 


BTI  Turbulence  Categories  (Lester  and  Burton, 


Cateoorv 

BTI 

Light 

60 

Moderate 

70 

Moderate/Severe 

90 

Severe 

100 

Extreme 

120 

1988) 
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(R  +  =  the  maximum  value  of  (R;  +  V,)  across 

the  model  domain. 

LTI  has  shown  promise  as  an  improved  method  for  LLT 
forecasts  in  the  Ft  Irwin  area,  and  is  the  basis  of  the 
system  developed  and  applied  in  the  following  sections. 
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3.  LLT  Algorithm  Development 

3 . 1  General 

Having  identified  the  BTI/LTI  as  a  basis  for  an  expert 
system  with  the  potential  for  meeting  both  scientific  and 
operational  needs  for  LLT  prediction  over  a  mesoscale  area, 
an  LLT  forecast  algorithm  for  Ft  Irwin  has  been  developed 
and  is  discussed  in  the  following  sections. 

The  forecast  system  is  designed  to  satisfy  the 
following  practical  requirements: 

(i)  User-friendly  to  compensate  for  differing  levels 
of  forecaster  computer  experience;  the  system  must  be  easily 
accessible  to  operators  of  all  experience  levels. 

(ii)  Executable  on  a  common  microcomputer  (e.g.,  IBM 
PC/AT  or  comparable  system) ;  the  system  actually  used 
employed  an  80286  microprocessor  and  80287  co-processor, 
with  a  40  megabyte  hard  disk  drive,  1  megabyte  of  random 
access  memory,  and  an  Extended  Graphics  Adaptor  (EGA)  color 
monitor. 

(iii)  Executable  in  nearly  real  time;  the  system  must 
provide  accurate  decision  assistance  in  a  timely  manner  to 
ensure  its  inclusion  in  the  forecast. 

(iv)  Minimal  keyboard  inputs,  because  forecaster  time 


is  often  at  a  premium. 
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(v)  Suitable  for  adaptation  to  other  areas  and  data 
bases. 

(vi)  Modular  structure  to  permit  improvements  to  the 
system  to  be  made  more  easily. 

In  order  to  meet  the  meteorological  requirements  listed 
above,  three  program  modules  are  required  to  process 
information  in  such  a  fashion  as  to  be  the  most  useful  for 
the  forecaster.  The  primary  forecast  decision  aid,  and  the 
part  which  required  the  major  development  effort,  is  the 
forecast  module.  Secondary  components  are  the  archive 
module .  which  stores  PIREPs  for  use  both  in  real  time  in  the 
forecast  module  and  in  future  research,  and  the  tutorial 
module,  which  instructs  the  forecaster  on  the  use  of  the 
program.  Figure  6  illustrates  the  program  structure  and  the 
relationships  of  the  modules  to  each  other.  Appendix  II 
presents  FORTRAN  code  for  all  elements  in  the  program. 

The  BTI  was  originally  developed  as  an  empirical, 
non-dimensional  index  derived  from  parameters  commonly 
available  from  aviation  observations  and  forecast  products. 
Consequently,  a  mix  of  units  (English  and  metric)  were  used. 
As  it  is  the  aim  of  this  program  to  allow  forecasters  to 
input  data  directly  from  available  resources,  the  mixed  unit 
convention  was  maintained  in  the  development  of  the  LTI,  and 
is  used  in  the  LTI  program  and  discussion  which  follows. 

All  necessary  unit  conversions  are  performed  in  the  computer 
code,  and  are  transparent  to  the  user. 
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3 . 2  Forecast  Module 

The  forecast  module  produces  the  actual  decision  aid. 

It  is  composed  of  six  FORTRAN  programs  and  one  C  program  for 
on-screen  graphics  (Figure  7) .  A  list  of  the  programs  and 
their  functions  is  presented  in  Table  6.  The  module  uses 
the  basic  procedure  developed  by  Lester  and  Burton  (1988)  to 
produce  an  initial  "raw"  Local  scale  Turbulence  Index  (LTI) 
at  1  km  intervals  across  the  NTC  complex.  The  "raw”  LTI  is 
modified  to  account  for  increased  "wake"  (lee  side) 
turbulence,  and  recent  PIREPs.  Contours  of  turbulence 
category  thresholds  are  displayed  superimposed  on  200  meter 
terrain  contours. 

3.2.1  Raw  Local  Turbulence  Index 

Roughness 

A  macroscale  (BTI)  roughness  of  59  (R  in  equation  1) 
was  calculated  from  the  difference  between  the  highest  and 
lowest  elevations  in  the  domain  of  the  model  (i.e.  Ft 
Irwin) .  The  LTI  roughness  values  at  each  grid  point  were 
then  computed  in  the  same  manner  for  a  1.0  km  square 
centered  on  the  grid  point.  All  roughness  values  were 
computed  using  100  meter  terrain  data  provided  by  the  US 
Aimy  Atmospheric  Sciences  Laboratory  (ASL) .  Contours  of 
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Figure  7  Forecast  Module  Architecture 
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Table  6  Forecast  Module  Architecture 


Program _ Function _ Files  Read  Files  Created 


FILEMAKR 

Keyboard  Input/ 
user  interface 

MAINDATA.DAT 

WINDMAKR 

Creates  WOCSS 
input  files 

MAINDATA.DAT 

WINDS . DAT 
TEMPS . DAT 

WOCSS 

Wind  model 

WINDS . DAT 
TEMPS . DAT 
TERRAIN.DAT 

S PD. DAT 

DIR.  DAT 

WINTER? 

Interpolates 
horizontal 
winds  to  61X61 
grid; 

calculates  w 

S PD. DAT 

DIR. DAT 

X. DAT 

Y. DAT 

WINDVEL.DAT 

WINDDIR.DAT 

W.DAT 

INDXMAKR 

Calculates  raw 
LTI 

MAINDATA.DAT 
WINDVEL.DAT 
WINDDIR.DAT 
RUFFILE . DAT 

RAWINDEX.DAT 

(WR.DAT) 

INDXMOD 

Modifies  LTI 
for  wake  and 
PIREPs . 

RAWINDEX.DAT 

W.DAT 

MAINDATA.DAT 
,  PIREP.DAT 
TERRTYPE . DAT 

NEWINDEX.DAT 

PRP.PLT 

LTIGRAPH 

Creates  plotter 
commands  for 
display 

’ntciooo.plt 

NEWINDEX.DAT 

PRP.PLT 

PLOT.DAT 

PLOT 

Graphics 

output 

PLOT.DAT 

Resident  data  files, 
j  Scratch  file. 

^  Created  in  ARCHIVE  module. 
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the  roughness  values  (based  on  1— kilometer  squares)  for  Ft 
Irwin  is  presented  in  Figure  8. 

Winds 

Wind  inputs  for  the  BTI  (V  in  Equation  1)  were 
initially  inferred  from  the  850  mb  surface  as  described 
below.  As  will  be  shown  in  a  later  section,  the  resulting 
BTI  values  were  unrealistically  high.  To  reduce  the 
magnitude  of  V,  and  hence  BTI,  wind  speed  at  5  m  was  taken 
from  a  representative  location  (Tower  2,  Figure  3).  See 
Section  4 . 1  for  a  description  of  the  meteorological  sensors 
used. 

LTI  (grid)  wind  inputs  were  calculated  using  the  Ludwig 
and  Endlich  (1988)  Winds  on  Critical  Streamline  Surfaces 
(WOCSS)  model.  The  WOCSS  model  uses  available  measures  of 
vertical  wind  and  temperature  profiles  to  define  surfaces 
within  which  air  flow  takes  place.  Critical  streamline 
concepts  are  used  to  define  these  surfaces,  which  can 
intersect  the  terrain  under  stable  conditions;  this  forces 
flow  around  topographical  obstacles.  In  this  application, 
WOCSS  provided  a  convenient  method  by  which  limited 
observations  could  be  objectively  assimilated  to  provide  a 
gridded  estimate  of  near-surf acc  winds.  Inputs  to  WOCSS 
include  terrain,  vertical  temperature  and  wind  profiles,  and 
surface  winds.  The  WOCSS  model  was  chosen  on  the  basis  of 
its  more  realistic  treatment  of  the  vertical  structure  of 


Figure  8  Contours  of  Ft  Irwin  Terrain  1-Kllometer  Roughne 

Values 

Roughness  contours  labeled  In  tens  of  feet.  Axes  are  UTM 
CIui)  coordinates.  See  text  for  roughness  calculation. 


40 


the  atmosphere  in  the  wind  interpolation,  and  because  it  was 
available  with  complete  documentation.  The  choice  of  wind 
interpolation  scheme  is  not  critical  to  the  mechanics  of  the 
LLT  forecast,  although  inaccuracies  here  will  lead  to  poor 
forecasts. 

While  the  WOCSS  model  is  capable  of  handling  a  61  x  61 
(i.e.,  1.0  km)  grid  point  domain,  the  resulting  executable 
file  is  1.9  megabytes,  too  large  for  most  microcomputers. 

Two  attempts  were  made  to  resolve  this  problem.  Initially 
the  nested  grid  capabilities  of  WOCSS  were  examined,  using  a 
course  grid  at  2.0  km  spacing  and  four  finer  grids,  each  at 
1.0  km  and  comprising  one  quarter  of  the  model  domain.  The 
four  resulting  1.0  km  wind  fields  were  then  merged  into  one. 
This  approach  proved  inadequate  due  to  inconsistencies  at 
the  boundaries  of  the  quarter  grids  and  excessive  time 
required  to  do  the  4  model  runs  (approximately  20  minutes) . 

The  second  approach  simply  calculated  winds  on  a  31  x 
31  (i.e.,  2.0  km)  grid  and  then  used  a  linear  interpolation 
to  obtain  the  1.0  km  values.  In  this  process,  the  initial 
winds  are  resolved  into  u  and  v  components,  interpolated, 
and  combined  into  direction  and  speed  values.  This  reduced 
the  time  to  complete  a  single  model  run  to  5  minutes. 

Towers  1-5  (Figure  3)  were  assigned  as  surface 
observation  Sites  1-5  (see  section  4.1)  in  the  WOCSS  input 
file  (Table  6) .  As  the  model  assigns  the  tower  locations  to 
the  nearest  grid  point,  it  became  necessary  to  alter  the 
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input  tower  elevations  to  match  the  grid.  This  prevented 
the  model  from  disregarding  tower  inputs  as  beneath  the 
terrain. 

For  operational  ease  and  due  to  sparse  data, 
observations  at  only  two  levels  were  used  to  specify  the 
vertical  structure  inputs  to  the  WOCSS  model;  surface  (i.e., 
Tower  2)  and  the  estimated  top  of  the  mixing  layer. 
Experience  has  shown  the  WOCSS  model  to  be  more  realistic 
when  the  top  of  the  domain  is  above  any  stable  layers  that 
will  affect  the  flow  (Ludwig,  1990) .  Since  terrain  heights 
vary  between  60  and  1780  m  above  mean  sea  level  (MSL) ,  a 
compromise  mixing  layer  height  of  2000  m  MSL  was  assumed, 
with  temperature  and  pressure  extrapolated  from  850  mb 
synoptic  charts  as  follows: 

Temperature  was  extrapolated  from  the  850  mb  level  to 
2000  m  using  the  surface  -  850  mb  lapse  rate.  The  lapse 
rate  was  calculated  using  temperatures  and  heights  from 
Tower  2  and  the  synoptic  850  mb  chart. 

Pressure  at  2000  m  was  estimated  from  the  synoptic 
850  mb  height  assuming  a  linear  10  mb  decrease  in  pressure 
per  100  m  increase  in  altitude.  This  assumption  is  a  valid 
approximation  near  850  mb  with  normal  850  mb  temperatures 
and  heights. 

As  2000  m  winds  are  not  normally  available,  they  must 
be  inferred  from  another  source.  Geostrophic  winds  computed 
from  surface  observations  are  inaccurate  due  to  the  complex 
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terrain  and  high  elevations.  Therefore,  winds  from  the 
nearest  synoptic  constant  pressure  surface  (850  mb)  were 
used  as  the  best  available  approximation  for  2000  m  winds 
for  input  to  the  WOCSS  model. 

Stability  and  Pressure  Tendency 

The  surface  -  850  mb  lapse  rate  calculated  for  the 
WOCSS  model  is  also  used  for  the  stability  input  to  the  BTI 
(S  in  Equation  1) .  Pressure  tendency  (T  in  Equation  1)  is 
taken  directly  from  Tower  2  output. 

3.2.2  Index  Modifications 

Once  a  "raw"  LTI  has  been  calculated  at  grid  points  via 
Equation  2,  the  program  adjusts  the  index  based  on  PIREPs 
and  anticipated  wake  enhancements  of  turbulence.  For 
aircrew  safety,  the  modifications  are  deliberately 
conservative,  i.e.,  the  raw  LTI  may  be  increased,  but  is 
never  reduced  by  the  index  modification  scheme. 

Wake  Enhancement 

This  modification  was  conceived  after  discussions  with 
Ft  Irwin  pilots  indicated  that  there  was  a  significant 
difference  in  turbulence  up-  and  downstream  of  mountain 
peaks  and  crests,  i.e.,  turbulence  was  more  common  on  the 
lee  slopes  of  the  mountains.  These  remarks  are  in  agreement 
with  experimental  and  theoretical  studies  of  airflow  over 
mountains  (e.g.,  Baines,  1987).  It  should  also  be  noted 
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that  neither  the  BTI/LTI  scheme  (Equations  1  and  2)  nor  the 
WOCSS  model  explicitly  deals  with  the  generation  of  wake 
turbulence . 

The  "lee"  slopes  of  individual  terrain  features  can  be 
identified  by  the  sign  of  the  terrain-induced  vertical 
velocities.  These  are  derived  from  the  dot  product  of  the 
interpolated  horizontal  wind  (V^)  and  the  topography  H(x,y) 
(Figure  9) ,  i.e. , 

w  =  ■  (6H/6xi  +  6H/6yj) 

where 

V„  =  horizontal  surface  wind  vector 
H  =  terrain  height  =  H(x,y) 

The  "lee”  sides  of  the  terrain  features  are  those  regions 
where  w  <  0  (Figure  9) . 

Assuming  that  wake  turbulence  is  more  likely  where 
magnitudes  of  vertical  velocities  are  large,  experiments 
were  conducted  to  determine  the  appropriate  w  threshold  at 
which  modify  the  LTI.  Thresholds  tested  were  w  <  0,  -0.5, 

-1  and  -1.5  ms  ’.  The  LTI  fields  generated  using  0  and  -0.5 
ms  ^  thresholds  were  unrealistically  large  compared  to  pilot 
observations,  i.e.,  they  resulted  in  the  selection  of  30-40% 
of  the  post  area  for  modification.  A  threshold  of  -1.5  ms’ 
was  more  realistic  and  is  in  agreement  with  an  observational 
study  by  Lester  (1974)  which  showed  that  a  vertical  velocity 
of  -1.5  ms  ’  was  the  threshold  for  "weak”  lee  waves  and 
their  associated  turbulence.  This  value  was  used  in  the 
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current  study.  An  example  is  shown  in  Figure  10. 

Once  the  areas  to  be  modified  were  identified,  the 
degree  of  modification  had  to  be  specified.  Consultation 
with  Burton  (1989)  indicated  that  forecast  turbulence  should 
increase  by  one  category  in  the  areas  identified  above.  The 
simplest  way  to  accomplish  this  is  to  increase  the  LTI  value 
by  20  in  these  areas.  This  has  the  desired  effect  of  an 
increase  of  one  turbulence  category  in  most  areas,  with  two 
exceptions:  (i)  areas  indicating  severe  turbulence  may  or 

may  not  increase  to  extreme,  a  moot  point  as  in  either  case 
aircraft  would  avoid  the  area,  and  (ii)  a  non-turbulent 
value  (<60)  could  be  upgraded  two  categories  to  moderate 
(>70) .  This  latter  case  is  unlikely,  since  the  horizontal 
wind  velocities  required  for  terrain  induced  vertical 
velocity  to  exceed  the  threshold  for  LTI  modification  are 
normally  large  enough  to  ensure  a  raw  LTI  greater  than  60. 
Figure  11  shows  the  program  logic  for  the  wake  modification. 

PIREP  Modification 

PIREPs  are  the  sole  connection  between  the  forecaster 
and  verification  of  his/her  turbulence  forecasts; 
consequently,  any  comprehensive  LLT  forecast  plan  must 
include  them.  The  method  for  incorporating  PIREPs  into  the 
LTI  is  described  in  this  section. 

Due  to  short  turbulence  time  scales,  current  forecast 
procedures  consider  PIREPs  useful  if  they  are  less  than  one 
hour  old  (AWS,  1988) .  On  this  basis,  and  after  the  raw  LTI 


Figure  10  Example  of  Wake  Modification  Criteria 


Polygonal  shaped  regions  are  areas  of  calculated  negative 
CdownwardD  vertical  winds  less  than  1.5  ms"*-.  Winds 
calculated  using  WOCSS  model  output  for  0000  GMT  1  March 


Figure  11  Wake  Modification  Program  Logic 


Raw  LTI  Is  increased  by  20  at  grid  points  with  negative 
vertical  velocities  greater  than  3  knots  In  magnitude. 
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has  been  modified  for  wake-enhanced  turbulence,  the  program 
searches  the  PIREP  data  base  (described  in  Section  3.3)  for 
reports  less  than  one  hour  old.  For  safety,  current  reports 
are  assigned  an  index  value  at  the  upper  end  of  their 
turbulence  category.  Therefore,  reports  of  "light," 
"moderate,"  and  "severe"  turbulence  would  be  assigned  LTI 
values  of  69,  89,  and  119  respectively  (Table  5) .  The  PIREP 
LTI  is  then  compared  to  the  LTI  at  the  nearest  grid  point. 

If  the  current  value  at  the  grid  point  is  greater  than  or 
equal  to  the  PIREP  value,  the  PIREP  is  ignored  and  the 
program  searches  for  more  PIREPs.  If,  on  the  other  hand, 
the  LTI  is  less  than  the  PIREP,  the  grid  point  is  assigned 
the  LTI  value  of  the  PIREP,  and  the  ratio  of  the  modified 
LTI  to  the  initial  LTI  is  calculated. 

After  modification  of  the  LTI  at  the  grid  point  nearest 
to  the  PIREP,  the  rest  of  the  LTI  field  is  considered  for 
modification  by  assuming  that  similar  terrain  in  other  areas 
will  induce  similar  turbulence  intensities.  Therefore  the 
LTI  field  is  modified  for  each  PIREP  in  all  areas  of  similar 
terrain. 

Although  terrain  can  be  described  by  a  number  of 
parameters,  such  as  altitude,  roughness,  slope,  orientation, 
etc.,  in  the  present  study,  slope  and  orientation  were 
chosen  as  the  most  important  terrain  characteristics  for  the 
production  of  turbulence  (Chapter  2) .  More  specifically, 
terrain  is  judged  to  be  similar  when  comparing  grid  points 
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if  the  direction  of  the  slope  is  within  30°  and  the 
magnitude  is  within  10  m  per  km.  These  data  have  been 
calculated  and  combined  into  a  four-digit  number  for  each 
grid  point  and  stored  in  a  terrain  characteristics  data 
file.  The  first  two  digits  represent  the  positive 
(up-slope)  direction  on  a  12-point  compass  (North  ±  15°  = 

1) .  The  final  two  digits  represent  the  magnitude  of  the 
terrain  slope  in  tens  of  meters  per  kilometer. 

During  program  execution,  the  terrain  characteristic  at 
the  PIREP  location  is  determined  first,  then  the  LTI  values 
at  grid  points  with  similar  terrain  (i.e.,  the  same  terrain 
characteristic  value)  are  multiplied  by  the  ratio  calculated 
earlier  (Figure  12) .  The  program  simply  searches  the 
terrain  characteristic  data  file  for  grid  points  with  values 
within  the  specified  ranges,  and  alters  the  LTI  at  those 
points. 


3.2.3  Display 

In  order  for  the  final  display  of  the  LTI  field  to  be 
useful,  it  must  be  unambiguous,  easily  read,  and  easily 
understood.  To  satisfy  these  requirements,  a  contour  map  of 
the  LTI  field  is  generated  on  the  computer  screen,  with 
turbulence  category  thresholds  (for  Category  I  aircraft) 
displayed  in  color.  200  m  terrain  contours  and  a  map  of 
Fort  Irwin  are  underlaid  in  the  background  for  reference. 


The  display  programs  were  provided  by  ASL  (Tabor, 
1988) ,  and  modified  for  use  with  the  LTI  program.  They 
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consist  of  one  FORTRAN  program  (LTIGRAPH) ,  and  one  C  program 
(PLOT) .  LTIGRAPH  is  a  contouring  program  that  creates  and 
stores  Hewlett-Packard  plotter  commands  in  a  data  file;  PLOT 
translates  these  to  screen  commands  for  the  EGA  color 
monitor.  While  PLOT  was  used  as  supplied,  LTIGRAPH  was 
modified  as  follows:  (i)  the  terrain  contours  were  fixed  at 
200  m,  (ii)  only  threshold  LTI  contours  were  displayed, 

(iii)  a  map  of  Ft  Irwin  boundaries  was  added,  (iv)  the 
resolution  was  increased  to  l.O  km,  and  (v)  provisions  were 
made  to  display  PIREP  locations  (denoted  as  "*”  in  the 
appropriate  color)  on  the  map. 

LTI  values  indicating  extreme  turbulence  cause  a 
warning  message  to  be  displayed  and  a  keyboard  input 
(carriage  return)  is  required  to  resume  LTIGRAPH  program 
execution.  This  section  of  the  program  was  used  as 
supplied,  although  the  threshold  for  this  message  is  easily 
changed  or  omitted. 

3 . 3  Archive  Module 


The  archive  module  creates  and  amends  the  PIREP  data 
file.  The  PIREP  data  file  is  a  direct  access  file  which 
stores  the  following  information  about  each  PIREP:  year. 
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day,  time,  location  (UTM) ,  turbulence  category,  and  aircraft 
category . 

The  forecast  module's  use  of  the  archive  module  as  a 
data  file  for  current  PIREPs  has  already  been  discussed. 

The  long-term  archival  of  these  data  is  an  equally  important 
function  of  this  module.  As  discussed  previously,  one  of 
the  primary  difficulties  in  developing  any  LLT  forecast 
program  is  lack  of  data  for  verification.  This  module  lays 
the  groundwork  for  future  studies  by  archiving  PIREPs;  it 
essentially  creates  the  data  base  needed  'or  future 
verification  or  modification  of  this  forecast  module,  and 
the  development  of  others. 

3 . 4  Tutorial  Module 

The  tutorial  module  was  not  developed,  although  it  is  a 
necessary  part  of  the  complete  system.  As  stated  before, 
this  module  should  support  new  forecaster  training  and  skill 
maintenance  for  experienced  forecasters.  To  that  end,  it 
should  contain  instructions  for  operation  (including 
examples) ,  and  sections  on  local  causes  of  LLT,  LTI 
background,  and  local  effects  (such  as  the  turbulence 
"climatology"  derived  by  Lester  and  Burton,  1988) . 
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4 .  Available  Data 

4 . 1  Tower  Data . 

In  order  to  create  and  test  forecast  methods, 
development  of  a  LLT  data  base  was  necessary.  Previous 
investigations  at  Ft  Irwin  have  been  hampered  by  the  lack  of 
data  for  verification.  As  noted  in  Chapter  1,  ASL  has 
recently  responded  to  the  problem  by  installing  instrumented 
towers  at  various  sites  throughout  the  reservation  (Marrs, 
1988) .  Lester  and  Burton  (1988)  have  also  initiated  a  PIREP 
collection  program  at  NTC.  These  data  are  now  available  for 
development,  initialization,  and  verification  of  the  LLT 
forecast  scheme  described  in  the  previous  section. 

Sensors  for  ASL's  Surface  Atmospheric  Measuring  System 
(SAMS)  were  mounted  on  five  5-meter  towers  in  various 
locations  throughout  the  Fort  Irwin  reservation  (Figure  3) 
chosen  by  ASL  on  the  basis  of  their  proximity  to  commonly 
used  operations  areas,  their  representation  of  "average" 
terrain  types,  and  their  accessibility  (Marrs,  1988) .  SAMS 
is  a  meteorological  data  collection  system  designed  and 
developed  for  ASL  by  the  New  Mexico  State  University 
Physical  Sciences  Laboratory.  Each  station  records 
pressure,  temperature,  relative  humidity,  wind  direction, 
wind  speed,  peak  wind  speed,  and  standard  deviation  of  the 
wind  direction.  During  the  collection  period,  all 
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information  was  transmitted  at  15-minute  intervals  to  the 
base  station  at  BYS. 

Printed  tower  data  for  the  period  of  interest  were 
manually  keyed  into  a  microcomputer  spreadsheet  program 
( Ouattro .  Borland  International)  for  use  in  the  current 
study.  Manual  transfer  was  necessary  as  prototype  system 
incompatibilities  precluded  electronic  transfer.  However, 
under  operational  conditions,  such  data  should  be  ingested 
automatically. 

4 . 2  PIREPs 

As  noted  previously,  Lester  and  Burton  (1988)  initiated 
the  collection  of  PIREPs  at  NTC.  Data  were  collected  over 
two  training  periods  between  27  February  and  12  May  1988, 
coinciding  with  the  SAMS  tower  study,  and  the  prevalence  of 
LLT  during  the  spring  transition  season. 

Pilot  reports  of  turbulence  were  collected  after  each 
mission  flown  during  the  two  cycles.  Because  of  their 
extensive  knowledge  of  the  Ft  Irwin  terrain  (Lester  and 
Burton,  1988) ,  only  permanent  party  pilots  were  asked  to 
participate  in  the  survey.  Each  pilot  was  requested  to  fill 
out  a  short  questionnaire  after  each  mission.  The  pilots 
were  asked  to  document  areas  of  turbulence  encountered  along 
their  flight  path,  including  date,  time,  type  of  turbulence 
encountered,  and  flight  level.  They  were  also  encouraged  to 
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report  areas  of  no  turbulence.  See  Appendix  III  for  a 
sample  PIREP  questionnaire. 

A  total  of  87  PIREP  forms  were  collected  and  screened 
for  usable  information.  Of  these,  approximately  20%  were 
unusable  for  a  variety  of  reasons  (e.g.,  missing  time,  date, 
route;  report  area  too  large,  etc.)  The  remaining  forms 
were  sorted  according  to  date  and  time.  Multiple  reports  on 
the  same  date  were  consolidated  onto  overlays  for  subsequent 
analysis.  Days  with  five  reports  or  more,  regardless  of 
time  of  day  or  turbulence  intensity,  were  selected  as  case 
days.  A  total  of  11  days  met  the  criteria  (Table  7) . 
Synoptic  conditions  for  these  case  days  have  been  documented 
by  Incerpi  (1989) .  A  total  of  19  sets  of  synoptic  data 
(0000  and/or  1200  GMT)  were  available  for  use  as  inputs  to 
the  program  (Table  8) . 
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Table  7 

PIREP  Case  Days 

Date 

SvnoDtic  Descriction 

29 

Feb  1988 

Pre-Frontal 

1 

Mar 

Post-Frontal 

5 

Mar 

Weak  Synoptic 

Gradient 

6 

Mar 

Weak  Synoptic 

Gradient 

22 

Apr 

Mountain  Wave 

25 

Apr 

Mountain  Wave 

26 

Apr 

Weak  Synoptic 

Gradient 

28 

Apr 

Pre-Frontal 

10 

May 

Post-Frontal 

11 

May 

Weak  Synoptic 

Gradient 

12 

May 

Weak  Synoptic 

Gradient 
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5.  Program  Execution  and  Results 

An  evaluation  of  the  forecast  model  is  presented 
below.  First,  a  detailed  example  of  a  model  run  is 
presented,  then  model  output  is  compared  to  observed  winds, 
LTI,  and  PIREPs  for  all  cases. 

5 . 1  Example  of  the  Program  Operation 

For  purposes  of  illustration  details  of  the  model 
steps  are  presented  below  for  a  single  case:  0000  hrs  GMT, 

1  March  1988.  It  should  be  noted  that  in  actual 
applications,  several  intermediate  steps  discussed  here  are 
transparent  to  the  user.  Only  the  inputs  and  final  displays 
(Appendix  IV)  are  presented  on  the  computer  screen. 

On  this  case  day,  the  reservation  was  under  the 
influence  of  a  high  pressure  system  to  the  southeast,  with 
post-frontal  conditions  and  moderate  southwesterly  winds. 
Table  9  shows  the  keyboard  inputs  for  the  case.  These  are 
controlled  by  the  FILEMAKR  program,  which  creates  the 
MAINDATA  data  file.  This  data  file  is  accessed  by 
subsequent  program  elements. 

After  the  FILEMAKR  keyboard  input  module,  the  WINDMAKR 
program  creates  two  data  files  for  use  in  the  WOCSS  code. 
Table  10  lists  the  data  in  these  files.  Once  these  files 
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Table  9  Keyboard  Inputs  for  1  Mar  1988,  0000  GMT 


Year  1988 

Month  3 

Date  1 

Time  (GMT)  0000 

Tower  2  Pressure  Tendency  (mb)  1.0 

BYS  (Tower  2)  Wind  Speed  (kt)  12 

850  mb  Wind  Speed  (kt)  20 

850  mb  Wind  Direction  210 

850  mb  Temperature  (°C)  10 

850  mb  Height  (m)  1495 


Tower  1 

Pressure  (mb)  919.5 

Temperature  (°F)  61 

Wind  Direction  236 

Wind  Speed  (kt) 


2 

3 

4 

5 

898.4 

892.2 

881.0 

934 

56 

58 

54 

64 

236 

202 

220 

227 

19 


12 


16 


13 


28 
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Table  10  Wind  Model  Inputs 

5  meter  data: 

Tower 


1 

2 

3 

4 

5 

UTM  East 

539.6 

521.5 

547.8 

539.5 

557.1 

UTM  North 

3898.6 

3895.8 

3927.5 

3920.2 

3901.7 

Pressure  (mb) 

919.5 

898.4 

892.2 

881.0 

934.6 

Temperature  (“F) 

61 

56 

58 

54 

64 

Wind  Direction 

236 

236 

202 

220 

227 

Wind  Speed  (kt) 

19 

12 

16 

13 

28 

Vertical  structure  (Tower  2  values  used  for  surface  input, 
2000  m  values  derived  as  shown  in  Section  3.2.1): 


Surface  2000  m 

(1026  m) 


UTM  East  521.5 

UTM  North  3895.8 

Pressure  (mb)  898.4 

Temperature  ( " C)  13.3 

Wind  Direction  ^  236 

Wind  Speed  (ms  ’)  6.2 


521.5 
3895.8 

799.5 
6.4 

210 

10.3 
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are  created,  the  WOCSS  model  is  executed.  Anemometer  height 
(10  m)  winds  are  saved  and  all  other  output  files  are 
deleted.  Figure  13  shows  model  output  winds  at  10  m.  Note 
the  lowest  height  wind  field  generated  by  the  WOCSS  model  is 
10  m;  these  are  used  as  an  approximation  to  the  tower  (5  m) 
winds  in  the  LTI  program. 

Once  the  horizontal  wind  field  is  generated  on  the  2  km 
grid,  the  WINTER?  module  uses  it  to  interpolate  1  km  winds, 
which  are  then  combined  with  the  resident  terrain  slope  data 
to  derive  the  w  field. 

The  INDXMAKR  routine  calculates  the  macroscale  BTI,  the 
(V+R) max  i^oiomalizing  factor  (Equation  2),  and  the  raw  LTI  at 
each  grid  point.  Figure  14  shows  the  raw  LTI  threshold 
values  and  200  m  terrain  contours.  The  raw  index  values 
range  from  35.5  (no  turbulence)  over  the  valleys  and  dry 
lakebeds,  to  111.0  (severe  turbulence)  over  the  mountain 
peaks . 

The  INDXMOD  program  adjusts  the  raw  LTI  for  wake  effects 
and  PIREPs,  and  creates  the  final  LTI  data  field.  Figure  15 
shows  the  wake-modified  LTI.  As  can  be  seen  by  comparison 
with  Figure  14,  the  modifications  coincided  with  areas  of 
large  negative  terrain  slope,  near  the  peaks  of  the  mountain 
ranges.  Note  the  maximum  LTI  value  has  increased  to  131.0 
(extreme  turbulence) ,  but  the  wake  modification  has  little 
effect  over  most  of  the  region. 
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Figure  13  WOCSS  Output  Tower  Winds 

10  meter  winds  generated  by  the  WOCSS  model  for 
0000  GMT  1  March  1988.  Arrow  at  lower  right  shows  scale  for 

10  ms"*^  wind. 
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Figure  14  Raw  CUnmodl f led^  LTI  200  m  terrain  contours  are 
shown  In  brown,  outline  of  NIC  In  black..  LTI  turbulence 
thresh-olds  are  outlined  according  to  the  legend  at  bottom 
right.  Axes  are  marked  In  2  km  Increments.  North  Is  at  the 
top  of  the  display. 
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Figure  15  Wake  Modified  LTI  As  for  Figure  14.  Arrows 
indicate  areas  of  modification. 
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In  order  to  illustrate  the  potential  impact  of  the 
PIREP  modification,  bogus  PIREPs  were  added  based  on  the 
turbulence  climatology  derived  by  Lester  and  Burton  (1988) . 
These  were  then  used  to  modify  the  LTI  with  the  results 
shown  in  Figure  16.  The  bogus  PIREPs  are  indicated  by 
asterisks,  and  are  concentrated  in  an  area  6  km  west  of  BYS 
(an  area  of  frequent  turbulence  encounters) .  Figure  16 
shows  the  extent  of  the  PIREP  modification,  with  many 
isolated,  individual  points  in  the  lake  beds  and  valleys 
being  modified  to  values  higher  than  those  in  Figure  15. 

The  modifications  are  more  continuous  over  the  mountain 
ranges.  The  maximum  LTI  value  has  increased  to  156.7.  Note 
th.’.s  was  for  demonstration  only;  bogus  reports  were  not  used 
fcr  the  actual  case  tests. 

5 . 2  Program  Results 

A  summary  of  observed  and  modelled  wind  speeds  and  LTI 
vr lues  is  presented  in  Appendix  V.  It  should  be  noted  that 
the  BTI  values  for  all  case  days  indicated  at  least 
moderate/ severe  turbulence  (Table  11) .  Consequently,  if  the 
BTI  were  the  only  diagnostic  tool  for  the  meteorologist,  the 
post  would  have  been  closed  to  flying  in  every  case.  Since 
BTI  has  been  used  widely  in  the  past,  this  result  is 
interpreted  as  typical  of  a  large  scale  turbulence  input. 

A  meaningful  evaluation  of  the  LTI  program  first 
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Table  11  BTI  Inputs  (R  =  59) 


DATE/TIME (GMT) 

T 

V 

S 

BTI 

29  FEB/1200 

5 

11 

27.1 

102 

01  MAR/0000 

10 

12 

29.9 

111 

01  MAR/1200 

12 

6 

25.7 

103 

05  MAR/1200 

4 

6 

24.9 

94 

06  MAR/0000 

18 

4 

34.3 

115 

06  MAR/1200 

13 

9 

24.8 

106 

07  MAR/0000 

13 

11 

33.4 

117 

22  APR/ 12 00 

6 

8 

30.0 

103 

23  APR/0000 

2 

17 

35.0 

113 

25  APR/ 12 00 

4 

9 

24.0 

96 

26  APR/0000 

15 

5 

35.5 

115 

26  APR/1200 

7 

2 

32.3 

100 

28  APR/1200 

7 

22 

29.9 

118 

10  MAy/1200 

1 

4 

28.6 

93 

11  MAY/0000 

18 

3 

36.6 

117 

11  MAY/1200 

2 

9 

31.1 

101 

12  MAY/0000 

17 

2 

40.6 

119 

12  MAY/1200 

3 

4 

30.8 

97 
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requires  an  examination  of  the  wind  model.  Observed  and 
modelled  winds  can  be  compared  at  the  tower  sites.  This  was 
accomplished  statistically,  and  results  are  presented  in 
Figure  17  for  all  19  cases.  The  modelled  wind  speeds  (V,,) 
are  related  to  the  observed  wind  speeds  (Vq)  as; 

=  1.09  Vp  +  1.98  ms’  (3) 

with  a  correlation  coefficient  of  0.70.  Equation  3  shows 
that  the  model  tends  to  overestimate  wind  speeds.  Figure  17 
also  shows  that  there  is  more  scatter  at  higher  wind 
speeds. ’ 

The  modelled  index  (LTI„)  at  each  of  the  towers  was 
also  compared  to  the  "observed”  index  (LTI^)  ,  i.e.,  the 
value  obtained  by  using  the  observed  tower  wind  and  1  km 
grid  roughness.  These  data  show  a  lower  correlation  than 
the  wind  speeds  (Figure  18) .  The  least-squares  fit  equation 
is: 

LTI^  =1.07  LTIq  -  4.61  (4) 

with  a  correlation  coefficient  of  .56.  It  is  noted  that  the 
model  shows  a  tendency  to  underestimate  at  low  LTI  values 
and  overestimate  at  higher  values,  reflecting  a  similar, 
though  lesser  magnitude,  tendency  in  the  WOCSS  model.  While 
any  discrepancy  is  unfortunate,  one  may  argue  that  the  LTI 


Since  these  tests  were  run,  the  WOCSS  code  has  been  modified 
so  that  winds  at  grid  points  nearest  the  input  sites  can  be 
constrained  to  remain  as  observed  (Ludwig,  1990) . 


Tower  Winds 

Observed  vs  WOCSS  Model 


5  10  15  20  25  30  35  40 

Observed  Wrd  (kts) 


17  Modelled  V3  Observed  Tower  Wl 
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errs  on  the  side  of  safety  for  the  higher  turbulence 
categories. 

Verifications  of  LTI  turbulence  estimates  for  all  cases 
were  made  by  comparing  them  with  PIREPs  as  described  in  the 
previous  chapter.  Tables  12  and  13  are  contingency  tables 
for  raw  and  wake-modified  LTI  vs  PIREPs  for  all  PIREPs 
within  +  1  hr  of  synoptic  times  (i.e.  case  times) .  Since 
pilots  reported  turbulence  areas  rather  than  precise 
locations,  LTI  values  were  determined  by  using  the  highest 
turbulence  category  underlying  the  reported  area  (see 
Appendix  III).  As  summarized  in  Table  14,  the  results  are 
nearly  identical  with  and  without  the  wake  modification. 

The  model  overestimated  turbulence  categories  52.2%  of  the 
time,  while  underestimating  only  4.3%.  Again,  the  errors 
are  on  the  side  of  safety. 

The  model  showed  better  correlations  with  PIREPs  for 
the  12  hr  period  following  the  model  run  (Tables  15  and  16) . 
Skill  scores  for  the  12  hr  forecast  were  .56  as  compared  to 
.29  for  the  nowcast  model. 

5.3  Discussion 

In  general,  the  LTI  overestimates  low  level  turbulence 
when  compared  with  PIREPs.  This  overestiroation  is  most 
likely  due  to  a  combination  of  overestimation  by  the  BTI,  by 
the  wind  interpolation  scheme  (already  documented) ,  and  the 
nature  of  the  verification  scheme  (i.e.,  PIREPs). 
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Table  13  Contingency  Table  of  Wake  Modified  LTI  vs  PIREPs 


OBSERVED 


FORECAST 

N 

L 

M 

S 

X 

TOTAL 

N 

9 

1 

10 

L 

1 

1 

2 

M 

1 

4 

5 

S 

2 

2 

4 

X 

2 

2 

TOTAL 

11 

10 

2 

23 
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Table  14  LTI  Nowcast  Statistics 


Raw 

Wake  Modified 

Number  of  Cases 

23 

23 

Total  Correct 

10 

10 

%  Correct 

43.5 

43.5 

%  Overestimate 

52.2 

52.2 

by  1  category 

30.4 

30.4 

by  2  categories 

17.4 

13.0 

by  3  categories 

4.3 

8.7 

%  Underestimate 

4.3 

4.3 

Sklll^ 

.29 

.29 

^Sklll  =  R  -  E 
T  -  E 


R  =  number  of  correct  forecasts 
T  =  total  number  of  forecasts 
E  =  number  of  forecasts  expected  to  be 
correct  based  on  chance 
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Table  15  Contingency  Table  of  Wake  Modified  LTI  vs  12-hr 

PIREPs 

OBSERVED 


FORECAST 

N 

B 

M 

S 

X 

TOTAL 

N 

33 

6 

B 

40 

L 

5 

14 

19 

M 

5 

5 

2 

12 

S 

2 

B 

2 

5 

X 

TOTAL 

45 

26 

5 

B 

76 

76 


Table  16  LTI  Forecast  Statistics 


Number  of  Cases  76 

Total  Correct  49 

%  Correct  64.5 

%  Overforecast  26.3 

by  1  category  15.8 

by  2  categories  7.9 

by  3  categories  2.6 

%  Underforecast  9.2 

by  1  category  7.9 

by  2  categories  1.3 

Skill  .56 
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An  attempt  to  compensate  for  the  BTI ' s  tendency  to 
overestimate  turbulence  categories  was  made  by  using 
near-surface  winds  instead  of  the  gradient  winds  which  were 
used  in  the  original  scheme.  Even  so,  the  BTI  values  for 
all  case  days  were  in  excess  of  the  threshold  of  90  required 
for  severe  turbulence.  This  may  indicate  that  the 
turbulence  category  thresholds,  which  were  originally 
derived  for  fixed-wing  light  aircraft  over  large-scale 
areas,  should  be  adjusted  for  helicopters  and/or  the  local 
terrain.  That  aircraft  were  allowed  to  fly  at  all  is 
indicative  of  the  lack  of  confidence  in  macro  scale  indices. 

The  WOCSS  wind  interpolation  scheme  showed  a  tendency 
to  overestimate  wind  speeds,  a  trait  which  may  be  attributed 
to  several  causes.  The  original  model  was  designed  for  use 
in  less  complex  terrain  than  is  encountered  at  NTC,  and 
might  not  handle  the  extremely  rough  terrain  encountered 
there.  Additionally,  the  number  of  towers  and  their 
locations  are  less  than  optimum;  they  may  provide 
unrepresentative  data.  More  towers  will  provide  more  input 
data,  and  presumably  a  better  model. 

Finally,  more  data  are  needed  for  verification.  Too 
few  PIREPs  were  available  for  a  statistically  meaningful 
test  of  the  LTI  algorithm.  Also,  in  many  cases,  the  areas 
marked  on  the  PIREP  questionnaires  were  obvious 
overestimates  of  the  actual  observation  area  (i.e.,  the 
circle  on  the  map  was  too  large) .  Consequently,  the  area 
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reported  may  cover  areas  of  high  LTI  which  were  not  actually 
overflown.  The  problem  of  the  subjective  nature  of  the 
PIREPs  has  already  been  addressed. 

The  model  performed  better  ’’forecasting"  LLT  in  the  12 
hours  following  the  model  run  than  it  did  in  describing  the 
current  situation.  As  the  preponderance  of  PIREPs  occurred 
between  1200  and  2359  GMT,  the  normal  diurnal  increase  in 
wind  speed  (mechanical  turbulence)  and  thermal  turbulence 
may  have  offset  the  tendency  toward  over-forecasts  in  the 
1200  GMT  model  run. 
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6.  Conclusions  and  Recommendations 

6.1  Conclusions 

The  Local  Turbulence  Index  (LTI)  was  applied  as  an 
objective  measure  of  turbulence  intensity  at  1  km  intervals 
across  NTC.  LTI  fields  showed  unambiguous  gradients  in 
turbulence  intensity  across  Ft  Irwin.  Typical  analyses 
which  showed  severe  turbulence  over  the  rugged  hills  in  the 
northern  portion  of  the  post  showed  lower  intensities  over 
the  broad  valleys  to  the  south. 

Two  modifications  to  the  original  formulations 
developed  by  Burton  (1964),  and  Lester  and  Burton  (1988) 
were  introduced.  In  the  first,  turbulence  areas  delineated 
by  certain  critical  LTI  values  were  enlarged  in  the  presence 
of  ''turbulent  wakes"  which  were  functions  of  the  local  wind 
velocity  and  the  slope  and  orientation  of  the  local  terrain. 
The  areal  impact  of  modification  was  small  because  the 
scales  of  the  wakes  are  small  compared  to  the  analysis  grid 
and/or  because  the  threshold  vertical  velocity  used  to 
define  the  wake  was  too  large.  Further  study  and 
verification  is  necessary. 

In  the  second  modification,  turbulence  reports  (PIREPs) 
taken  within  one  hour  of  the  LTI  analysis  were  assigned  an 
LTI  value  representative  of  the  reported  intensity. 
Subsequently,  the  slope  and  orientation  of  the  terrain  was 


used  to  identify  "similar”  locations  across  NTC;  the  LTI  at 
the  latter  points  were  modified  to  the  same  degree  as  the 
initial  point.  The  areal  impact  of  this  modification  was 
significant.  However,  there  was  some  question  as  to  the 
reality  of  the  pattern  (i.e.,  many  small,  closed  "turbulent" 
regions  appeared  across  the  domain) .  It  is  suggested  that 
the  addition  of  terrain  height  and  local  wind  velocity  to 
terrain  slope  and  orientation  would  be  more  useful  in  the 
search  for  "similar"  locations.  Further  verification  with 
more  and  better  PIREPs  is  needed. 

The  algorithm  developed  here  is  a  useful 
nowcast/ forecast  system  for  LLT  not  associated  with  moist 
convection.  It  meets  the  definitions  of  an  "expert  system" 
presented  in  Chapter  2.  It  yields  a  consistent  product, 
eliminating  forecaster  bias  and  bringing  all  forecasters  to 
a  similar  experience  level.  The  forecaster  need  not  be 
familiar  with  the  local  terrain  or  conditions  in  order  to 
provide  reasonably  accurate  forecasts  of  LLT.  The  use  of  a 
numerical  index  eliminates  the  subjectivity  associated  with 
descriptive  categories  and  allows  the  effective  use  of 
statistics  in  the  verification  and  improvement  of  the 
program. 

The  system  is  operationally  oriented  and  is  designed 
for  use  on  the  type  of  microcomputer  commonly  found  in 
forecast  offices.  It  requires  minimal  operator  inputs,  with 
the  exception  of  the  tower  data.  The  time  required  to 
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execute  the  model  is  on  the  order  of  10  minutes,  with  most 
of  that  time  being  spent  on  wind  model  computations  (i.e., 
forecaster  "hands  off”  time) . 

The  LTI ' s  usefulness  as  an  operational  LLT  forecast 
tool  is  the  result  of  several  characteristics.  The  index 
successfully  resolves  the  conflict  in  scale  between  point 
(mesoscale)  warnings  issued  for  the  area  and  the  microscale 
environment  at  Ft  Irwin.  Whereas  the  BTI  values  for  all 
case  days  were  in  excess  of  the  threshold  of  90  required  for 
severe  turbulence,  the  LTI  indicated  that  the  severe 
turbulence  was  usually  confined  to  small  areas  over  the 
mountain  tops.  The  LTI  showed,  and  PIREPs  verified,  the 
absence  of  severe  turbulence  over  most  of  the  reservation, 
emphasizing  the  inadequacy  of  a  point  warning  approach. 

The  graphic  display  of  the  LTI  field  is  the  ideal  way 
to  present  these  data  for  interpretation.  The  usefulness  of 
the  display  is  limited,  however,  by  the  lack  of  a  paper  copy 
backup  to  facilitate  pilot  briefings  and  provide  a  pemmanent 
record  of  the  program  output. 

The  comparison  of  the  model  as  a  forecast  and  nowcast 
aid  with  observations  illustrates  its  utility,  as  well  as 
its  tendency  to  overestimate  turbulence  intensity.  This 
overestimation  could  result  from  a  number  of  sources, 
including  the  wind  model,  the  scaling  procedure 
(Equation  2) ,  modification  procedures,  and  the  LTI 
turbulence  category  threshold  values. 
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The  archive  module,  if  properly  maintained,  will  prove 
extremely  valuable  for  future  investigations  of  LLT  at  Ft 
Irwin  by  providing  a  comprehensive  data  base  of  turbulence 
obseirvations . 

6.2  Recommendations 

Recommendations  for  further  study  and  development  of 
the  program  fall  into  three  categories:  (i)  improvements  to 
execution  of  the  current  program,  (ii)  improvements  to  the 
forecast  methodology,  and  (iii)  improved  applications.  Of 
these,  the  improvements  to  program  execution  can  be  realized 
most  rapidly. 

Short  term  improvements  will  reduce  the  workload  on  the 
forecaster  and  improve  the  display  of  the  final  index.  For 
example,  a  subroutine  to  read  the  tower  data  directly  from 
the  base  station  computer  is  required,  and  to  this  end  the 
section  is  isolated  in  the  code.  As  stated  in  the  previous 
section,  a  paper  copy  of  the  final  output  field  will  be 
beneficial.  The  illustrations  of  the  screen  display 
presented  here  were  plotted  on  ?  Hewlett-Packard  HP  7475A 
plotter;  however  plot  times  were  on  the  order  of  30-45 
minutes.  An  operational  forecaster  would  not  have  that  much 
time.  The  alternative  is  to  use  a  dot  matrix  printer,  but 
this  poses  the  problem  of  superimposing  many  different 
contours  in  black  and  white.  A  possible  solution  would  be 
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to  contour  the  terrain  and  only  one  critical  LTI  threshold, 
indicating  fly/no  fly  areas. 

More  data  are  needed  for  testing  and  modification  of 
this  and  other  LLT  forecast  systems,  both  at  Ft  Irwin  and  in 
other  areas.  To  this  end,  a  major  effort  should  be  made  to 
gather  quantitative  turbulence  data,  such  as  that  provided 
by  aircraft-borne  accelerometers,  although  the  PIREP  data 
base  should  not  be  neglected.  Improvements  to  the  forecast 
methodology  will  also  have  to  wait  for  a  larger  turbulence 
data  base. 

Clearly,  the  accuracy  of  the  wind  interpolation  scheme 
plays  a  major  role  in  determining  the  accuracy  of  the  LTI. 

As  has  been  discussed,  numerous  wind  models  are  available. 
Operationally,  a  compromise  must  be  made  between  speed 
(i.e.,  model  simplicity)  and  accuracy.  A  simpler  model 
might  be  able  to  interpolate  winds  to  the  full  61  x  61  grid, 
which  may  be  more  accurate  in  the  long  run  as  it  would 
eliminate  linear  interpolation  on  the  31  x  31  grid.  The 
assumptions  made  in  the  vertical  structure  inputs  to  the 
WOCSS  model  might  also  be  avoided  with  a  simpler  model. 

A  review  of  the  LTI  wake  and  PIREP  modification 
routines  in  the  light  of  a  larger  turbulence  data  bases 
holds  promise  for  great  improvements.  The  wake  and  PIREP 
modifications  both  involved  arbitrary  decisions  based  on 
research  and  past  experience  which  should  be  examined  in  the 
light  of  future  data.  The  wake  modification  may  be  more 


appropriately  applied  at  a  different  vertical  velocity 
threshold,  or  perhaps  in  conjunction  with  other  parameters 
such  as  distance  from  major  ridgelines,  stability,  etc. 

Also  the  actual  modification  might  be  better  applied  as  a 
percentage  change  in  the  index  (as  in  the  PIREP  case) . 

Another  weak  link  in  the  PIREP  modification  scheme  is 
in  the  determination  of  where  else  to  apply  the 
modification.  Different  terrain  param.eters  should  be 
examined  for  use  as  selection  criteria.  Also,  other 
modifications  should  be  considered,  such  as  the  inclusion  of 
a  low  level  wind  shear  parameterization. 

Other  aspects  of  this  program  that  should  be  examined 
further  with  respect  to  its  application  include  validity  of 
turbulence  category  thresholds  for  the  LTI,  flight  levels 
for  which  the  index  is  valid,  and  "transportability"  to 
other  areas.  The  question  of  appropriate  terrain  scale,  and 
hence  model  grid  spacing,  should  also  be  addressed.  The 
answer  will  likely  involve  a  compromise  between  optimum 
resolution  and  cost  (in  terms  of  both  time  and  computer 
power) . 

The  LLT  nowcast/ forecast  program  presented  is  a  useful 
tool  which  can  be  used  in  its  current  configuration,  and  can 
be  improved  with  the  development  of  a  larger  PIREP  data 
base.  It  is  the  author's  hope  that  the  program  will  serve 
as  a  foundation  for  future  studies. 
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APPENDIX  I: 


Turbulence-Prone  Areas  at  Ft  Irwin 


The  map.  redrawn  from  Lester  and  Burton  C1988), 
outlines  turbulence  prone  areas  of  Ft  Irwin.  It  is  part  of 
the  turbulence  "climatology"  derived  from  Interviews  with 
permanent  party  pilots  stationed  at  Ft  Irwin.  For  more 
detailed  descriptions  of  the  areas  see  Lester  and  Burton 
C1988D. 


nnnn 


94 


PROGRAM  ARCHIV 


THIS  PROGRAM  LOADS  PILOT  REPORTS  INTO  THE  PIREP.DAT 
DATA  FILE. 


CHARACTER  INTENS,  ANS 
INTEGER  CDAY.ACCAT 
REAL  MONTH 
LOGICAL  THERE 
CALL  CLEAR 

WRITEC *.* D ' Please  answer  the  following  questions  about 

» 

WRITEC * . * D ' the  turbulence  incident:  ' 

WRITEC*. O 
WRITEC* , 

200  WRITEC*. *)’What  was  the  YEAR  CGMT.  four  digits)?  ' 
READC*  .  *)YE?»R 

WRITEC* , *) 

IFCYEAR.lt. 1988. )THEN 
WRITEC* ,3) 

GO  TO  200 
END  IF 

201  WRITECS *)'What  was  the  MONTH  CGMT,  two  digits)?  ' 
READC*. ‘)MONTH 

WRITEC*. •) 

IFCCMONTH.GT. 12.). OR. CMONTH.lt. 1. ))THEN 
WRITECS3) 

GO  TO  201 
END  IF 

202  WRITEC* , *)' What  was  the  DATE  CGMT,  two  digits)?  ' 
READC*. •)DAY 

WRITEC* , •) 

IFCCDAY.lt. 1. ). OR. CDAY.GT. 31. ))THEN 
WRITEC*, 3) 

GO  TO  202 
END  IF 

CALL  DATECDAY, MONTH. YEAR, CDAY) 

203  WRITEC* ,*)' What  was  the  TIME  CGMT,  four  digits)?  ' 
READC* , *)TIME 

WRITEC*. *) 

IFCCTIME.LT.O. ).0R.CTIME.GT.2359. ))THEN 
WRITEC*. 3) 

GO  TO  203 
END  IF 

204  WRITEC*, •)' Enter  the  UTM  COORDINATES  in  km.  ' 

WRITEC* , *) 

WRITEC*,*)'  NORTH;  ' 

READ  C*. •)  y 
WRITEC* , *) 

WRITEC*,*)'  EAST:  ’ 

READ  C* , *)  X 


an  <  o  n  o 
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WRITEC‘, 

IFCCY.lt. 3887 . D . OR . CY . GT . 3947 . ) . OR. CX . LT . 506 . ID . OR . 

+  CX.GT.566. IDDTHEN 

WRITEC*,3D 
GO  TO  204 
END  IF 

205  WRITEC* , ‘D'What  was  the  reported  TURBULENCE  INTENSITY? 


WRITEC* , *D'  L  /  M  /  S  /  X  ' 

READ  C*.2D  INTENS 
WRITEC*. *D 

IFCC INTENS. NE. ’ L' D . AND . C INTENS. NE . ' 1' DDTHEN 
IFCCINTENS.NE. 'M' D. AND. C INTENS. NE. ' m' DDTHEN 

IFCCINTENS.NE. 'S' D. AND. C INTENS. NE. ' s' DDTHEN 
IFCC INTENS. NE. ' X' D . AND . C INTENS . NE . 'x' DDTHEN 
WRITEC*. 3D 
GO  TO  205 
END  IF 
END  IF 
END  IF 
END  IF 

206  WRITEC*, ‘D'What  was  the  AIRCRAFT  CATEGORY?  Cl, 2, 3D  ' 
READ  C* , *D  ACCAT 
WRITEC*, *D 

IFCCACCAT . LT . ID . OR. CACCAT . GT . 3DDTHEN 
WRITEC* , 3D 
GO  TO  206 
END  IF 

CALL  CLEAR 
ERRROR  CHECK 

WRITEC *,* D ' You  have  entered  the  following  values;' 
WRITEC *,4D  YEAR 
WRITEC *,5D  MONTH 
WRITEC*, 6D  DAY 
WRITEC*, 7D  TIME 
WRITEC*, 8D  Y 
WRITEC*, 9D  X 
WRITEC *,20D  INTENS 
WRITEC*, 21D  ACCAT 
WRITEC* , *D 

WRITEC*, *D  'Press  "RETURN”  to  continue,  "U"  to  update 
alues  ' 

READC* ,2DANS 

IFCCANS.EQ. 'u' D.OR.CANS.EQ. 'U' DDTHEN 
CALL  CLEAR 
GO  TO  200 
END  IF 

AAA***AA*AAA*AAAA****AAAA*A*A**A* 

CHECK  FOR  PI REP  FILE 


ooooo  o o o  •  ooo  nno  noo 
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INQUIRECFILE=’PIREP.DAT' ,  EXIST  >=  THERE> 

IFC.NOT.THERED  CALL  MPRP 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

OPEN  PI REP  FILE  AND  READ  NUMBER  OF  RECORDS 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

0PENC10,FILE=' PIREP.DAT' , ACCESS^' DIRECT' , STATUS=' OLD' , 
+  FORM=' FORMATTED' ,  RECL=50) 

ICOUNT-0 
DO  10,  1=1,999 

READCIO, 101,REC=I,END=11)  CHECK 
IC0UNT*IC0UNT+1 

10  CONTINUE 

11  IC=IC0UNT+1 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

WRITE  TO  FILE 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

WRITECIO, 100, REC=IC)YEAR,CDAY, TIME, X,Y, INTENS, ACCAT 
CLOSECIOD 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CHECK  FOR  MORE  REPORTS 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CALL  CLEAR 

WRITEC ‘ * D ' Do  you  wish  to  record  another  PIREP?  CY/N) 


READ  C*,2)  REPLY 
IFCCREPLY.EQ. 'N' D.OR.CREPLY.EQ. 
GO  TO  200 

2  FORMATCAID 

3  FORMATCIX, 

4  FORMATCIX, 

5  FORMATCIX, 

6  FORMATCIX, 

7  FORMATCIX, 

B  FORMATCIX, 

9  FORMATCIX, 

20  FORMATCIX, 

21  FORMATCIX, 

100  FORMATCIX, F6 

101  F0RMATCF6.0D 
9999  END 


n'))GO  TO  9999 


INPUT  BEYOND  LIMITS',/) 
Year;  ',F5.0) 
Month: 

Day; 

Time : 

Location: 


,F3.0) 

,F3.0) 


Intensity: 
A/C  Category; 
0,2X, 13, F6. 


' ,F5.0) 

' ,F7. 1, ' 

' ,F6. 1, ' 

'  ,A1) 

'  ,11) 

0,2X,2CF11 


north’  ) 
east '  ) 


1,2X), A1,2X, II) 


SUBROUTINE  CLEAR 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  TO  CLEAR  THE  SCREEN 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


DO  10,  1=1,30 
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WRITEC* . 

10  CONTINUE 
RETURN 
END 
C 
C 

C  - 

SUBROUTINE  DATECDAY, MONTH, YEAR, CDAYD 
C  - 

C 

Q  *A**A*ll**«**«*«Jk**A***lk**A******* 

C  THIS  S/R  CALCULATES  JULIAN  CCALENDAR:)  DATES  FROM  DAY, 
MONTH, 

C  YEAR  VALUES  AND  WAS  DEVELOPED  BY  F.D.  BARLOW  CBARLOW, 
1980, 

C  'SOLAR  ENERGY'  25,  P479).  MODIFICATIONS  BY  N.L. 

SMITH,  1988. 

REAL  DAY, MONTH, YEAR 
INTEGER  CDAY 
YRES  =  AM0DCYEAR,4.0D 
CRES  =  AMODCYEAR, 100.0) 

FACTOR  *  32.8 

IF  C YRES. EQ. 0.0. AND. CRES. NE.0.0)FACT0R  =  31.8 

IF  CM0NTH.LE.2.0)FACT0R  *  30.6 

CDAY  »  INTC30.6*M0NTH+DAY-FACT0R+.5) 

RETURN 

END 

C 

C 

C  - 

SUBROUTINE  MPRP 

C  - 

C  INITIALIZES  PIREP  FILE 

C  - 

CHARACTER  INTENS 

0PENC10,FILE=’PIREP.DAT' , ACCESS^' DIRECT' , STATUS=' NEW' . 
+  FORM*' FORMATTED' ,  RECL=50) 

YEAR“9999. 

IDAY-999 
TIME=999.9 
X*99999.9 
Y»999999.9 
INTENS*' Z' 

ICCAT-9 

WRITECIO, 100,REC=1)YEAR, IDAY , T IME, X , Y , INTENS, ICCAT 
100  F0RMATC1X,F6.0,2X, I3,F6.0,2X,2CF11. 1,2X), A1,2X, II) 
CLOSECIO) 

RETURN 

END 


n  o  o 
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200 


201 


202 


203 


?  ' 


PROGRAM  FILEMAKR 

********************************* 

CREATES  DATA  FILE  FOR  USE  WITH  LLT  FORECAST  AID 
**>1****************************** 


CHARACTER  ANS 

INTEGER  CDAY.TIME 

REAL  DAY , MONTH, YEAR. LAPSE 

REAL  TPRESC5) , TTEMPC5) , TD IRC5) , TVELC5) , M2FT 
CALL  CLEAR 

WRITEC* , What  is  the  YEAR  CGMT,  four  dlgltsD?  ’ 
READC* , *)YEAR 
WRITEC*. *3 

IFCYEAR.lt . 1988 . 3THEN 
WRITEC*,33 
GO  TO  200 
END  IF 

WRITEC* , *3' What  is  the  MONTH  CGMT,  two  dlglts3?  ' 
READC*. ‘3M0NTH 
WRITEC* , *3 

IFCCMONTH.GT. 12. 3.0R. CMONTH. LT. 1 . 33THEN 
WRITEC‘.33 
GO  TO  201 
END  IF 

WRITEC* . *3'What  Is  the  DATE  CGMT,  two  dlglts3?  ' 

READC* . •3DAY 
WRITEC*. ‘3 

IFCCDAY.lt. 1. 3.0R.CDAY.GT.3.. 33THEN 
WRITEC* ,33 
GO  TO  202 
END  IF 

CALL  DATECDAY, MONTH, YEAR, CDAY3 

WRITEC* , *3' What  is  the  TIME  CGMT,  four  dlglts3?  ' 
READC*. *3TIME 
WRITEC* , *3 

IFCCTIME.LT.O. 3.0R.CTIME.GT.2359. 33THEN 
WRITEC*. 33 
GO  TO  203 
END  IF 

WRITEC  *,*  3 '  What  is  the  three-hour  pressure  change  Cnib3 


READ  C* , *3  APP 
WRITEC* , *3 

WRITEC *,* 3 ' What  Is  the  surface  wind  speed  In  knots?  ' 
READ  C*,*3  BYSWIND 
WRITEC* , *3 
CALL  CLEAR 

WRITEC *,* 3 ' You  have  entered  the  following  values:  ' 

WRITEC*, *3 

WRITEC*, 43  YEAR 

WRITEC *,53  MONTH 

WRITEC* ,63  DAY 


o  o  o 
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WRITEC*.7)  TIME 
WRITEC*.8)  APR 
WRITEC*.9D  BYSWIND 
WRITEC*, 

WRITEC*,*)  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READCS2)ANS 

IFCCANS.EQ. 'u' D.OR.CANS.EQ. 'U' )DTHEN 
CALL  CLEAR 
GO  TO  200 
END  IF 

CALL  CLEAR 

204  WRITEC * . * D ' For  850  mb,  what  Is  the' 

WRITEC* , ‘D 

WRITEC* I  * Wind  speed  Cknots)  ?  ' 

READ  UAVEL 

WRITEC* , 

WRITEC  *  >  *  D '  Wind  direction  Cdeg!)  ?  ' 

READ  C* ,  *:>  UADIR 
WRITECS 

WRITEC D ' Temperature  Cdeg  C)  ?  ’ 

READ  C*.*D  UATEMP 
WRITEC*  ,  *:) 

WRITEC- Height  CmetersD  ?  ' 

READ  UAHT 

WRITEC* , *3 
CALL  CLEAR 

WRITEC *»* D' You  have  entered  the  following  values  for 
850  mb:  ' 

WRITEC* . *) 

WRITEC *,10)  UAVEL 
WRITEC*, 11)  UADIR 
WRITEC*, 12)  UATEMP 
WRITEC*, 13)  UAHT 
WRITEC* , *) 

WRITEC*,*)  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READC* ,2)ANS 

IFCCANS.EQ. 'u' ). OR. CANS. EQ. 'U' ))THEN 
CALL  CLEAR 
GO  TO  204 
END  IF 

CALL  CLEAR 


BEGIN  TOWER  INPUT 


CALL  CLEAR 

205  WRITEC* , *)'What  is  the  PRESSURE  Cm  mb)  at  Tower  1  ?  ' 
READ  C* , *)  TPRESCl) 

WRITEC* , *) 
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WRITEC* . *D' What  la  the  TEMPERATURE  Cm  “F)  at  Tower  1 

?  ' 

READ  C*,*)  TTEMPCID 
WRITEC*, 

WRITEC*, ‘D'What  la  the  wind  DIRECTION  Cm  "D  at  Tower 

1  ?  ' 

READ  C*,*D  TDIRCID 
WRITEC*. 

WRITEC*. ‘D'What  la  the  wind  VELOCITY  at  Tower  1  ?  ' 
READ  C* . *)  TVELCID 
CALL  CLEAR 

WRITEC * » * ) ' You  have  entered  the  following  values  for 
Tower  1:  ' 

WRITEC* . O 

WRITEC‘,14D  TPRESCi:) 

WRITECS15D  TTEMPCID 
WRITEC*.11D  TDIRCID 
WRITECNIOD  TVELCl) 

WRITEC*. 

WRITEC*,*D  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READC* .2DANS 

IFCCANS.EQ. ' u' ).0R.CANS,EQ. 'U' ))THEN 
CALL  CLEAR 
GO  TO  205 
END  IF 

CALL  CLEAR 

206  WRITEC* , *)'What  Is  the  PRESSURE  Cm  mb)  at  Tower  2  ?  ' 
READ  C‘.*)  TPRESC2) 

WRITEC* . *) 

WRITEC* , *)' What  is  the  TEMPERATURE  Cm  “F)  at  Tower  2 

?  ' 

READ  C* , *)  TTEMPC2) 

WRITEC* . *) 

WRITEC *,*)' What  Is  the  wind  DIRECTION  Cln  “)  at  Tower 

2  ?  ’ 

READ  C*, *)  TDIRC2) 

WRITEC  * , * ) 

WRITEC*! •)' What  la  the  wind  VELOCITY  at  Tower  2  ?  ' 
READ  C*,*)  TVELC2) 

CALL  CLEAR 

WRITEC *,*)' You  have  entered  the  following  values  for 
Tower  2;  ' 

WRITEC*. *) 

WRITEC*, 14)  TPRESC2) 

WRITEC*, 15)  TTEMPC2) 

WRITEC*, 11)  TDIRC2) 

WRITEC*, 10)  TVELC2) 

WRITEC* . *) 

WRITEC*,*)  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 
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READC‘,2)ANS 

IFCCANS.EQ. 'u' D.OR.CANS.EQ. 'U' DDTHEN 
CALL  CLEAR 
GO  TO  206 
END  IF 

CALL  CLEAR 

207  WRITEC* ,  O'What  is  the  PRESSURE  Cm  mb)  at  Tower  3  ?  ' 
READ  C‘.*D  TPRESC3D 

WRITEC*, 

WRITEC* , What  is  the  TEMPERATURE  Cln  “F)  at  Tower  3 

?  ' 

READ  C*,*D  TTEMPC3D 
WRITEC*, 

WRITEC*. O'What  Is  the  wind  DIRECTION  Cln  at  Tower 

3  ?  ' 

READ  C*, ‘D  TDIRC3D 
WRITEC*>  *3 

WRITEC*, *3'What  is  the  wind  VELOCITY  at  Tower  3  ?  ' 
READ  C*,*3  TVELC33 
CALL  CLEAR 

WRITEC* ,* 3' You  have  entered  the  following  values  for 
Tower  3;  ' 

WRITEC*, *3 

WRITEC*, 143  TPRESC33 
WRITEC*, 153  TTEMPC33 
WRITEC*, 113  TDIRC33 
WRITEC*, 103  TVELC33 
WRITEC*, *3 

WRITEC*, *3  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READC* ,23ANS 

IFCCANS.EQ. 'u' 3.0R. CANS.EQ. 'U'33THEN 
CALL  CLEAR 
GO  TO  207 
END  IF 

CALL  CLEAR 

208  WRITEC* , *3' What  is  the  PRESSURE  Cln  mb3  at  Tower  4  ?  ' 
READ  C*,*3  TPRESC43 

WRITEC* , *3 

WRITEC* ,  *3' What  is  the  TEMPERATURE  Cln  '’F3  at  Tower  4 

?  ' 

READ  C*,*3  TTEMPC43 
WRITEC*, *3 

WRITEC* , *3' What  Is  the  wind  DIRECTION  Cln  “3  at  Tower 

4  ?  ' 

READ  C*  ,  *3  TDIRC'O 
WRITEC*, *3 

WRITEC* , *3' What  is  the  wind  VELOCITY  at  Tower  4  ?  ' 
READ  C*,*3  TVELC43 
CALL  CLEAR 


o  o  n 
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WRITEC  *  »  *  !) '  You  have  entered  the  following  values  for 
Tower  A:  ' 

WRITEC*, 

WRITEC*, lAD  TPRESCAD 
WRITEC*, 15)  TTEMPCA) 

WRITEC*, 11)  TDIRCA) 

WRITEC*, 10)  TVELCA) 

WRITEC*, *) 

WRITEC*,*)  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READC* ,2)ANS 

IFCCANS.EQ. 'u' ).0R.CANS.EQ. 'U' ))THEN 
CALL  CLEAR 
GO  TO  208 
END  IF 

CALL  CLEAR 

209  WRITEC* , *)'What  Is  the  PRESSURE  Cm  mb)  at  Tower  5  ?  ' 
READ  C*,*)  TPRESC5) 

WRITEC*, *) 

WRITEC*, *)' What  Is  the  TEMPERATURE  Cm  “F)  at  Tower  5 

?  ’ 

READ  C*,*)  TTEMPC5) 

WRITEC* , *) 

WRITEC* ,*)' What  Is  the  wind  DIRECTION  Cln  °)  at  Tower 

5  ?  ' 

READ  C*,‘)  TDIRC5) 

WRITEC* , *) 

WRITEC*, *)'What  is  the  wind  VELOCITY  at  Tower  5  ?  ' 
READ  C*,*)  TVELC5) 

CALL  CLEAR 

WRITEC *,*)' You  have  entered  the  following  values  for 
Tower  5:  ' 

WRITEC* , *) 

WRITEC *,1A)  TPRESC5) 

WRITEC *,15)  TTEMPC5) 

WRITEC*, 11)  TDIRC5) 

WRITEC*, 10)  TVELC5) 

WRITEC  * , * ) 

WRITEC*’,*)  'Press  "RETURN"  to  continue,  "U"  to  update 
values  ' 

READC* ,2)ANS 

IFCCANS.EQ. 'u' ).0R.CANS.EQ. 'U' ))THEN 
CALL  CLEAR 
GO  TO  209 
END  IF 

CALL  CLEAR 

*AA*A***ftAft«ft*ft*ll*********  ******* 

CALCULATE  LAPSE  RATE  CDEG  C  PER  1000  FT) 

********************************* 

M2FT  =  3.2808399 
DT  =  TTEMPC2)  -  UATEMP 


nan  non 
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DZ  “  CUAHT  -  1026D  *  M2FT 
LAPSE  =  ABSCDT/DZ  *  1000) 

APP  -  10  ‘  ABSCAPP) 

A******************************** 

OPEN  FILE  AND  WRITE  DATA 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

OPEN  CUNIT  =  10. 

+  FILE  =  'MAINDATA.DAT', 

+  ACCESS  “  'DIRECT', 

+  STATUS  =  'NEW', 

+  FORM  “  'FORMATTED', 

+  RECL  =  10) 

WRITEC10.2000,REC=1)YEAR 
WRITECIO, 1000,REC=2)CDAY 
WRITECIO, 1000, REC=3)TIME 
WRITECIO, 2000, REC=4)APP 
WRI TEC 10 , 2000 . REC=5)BYSW I ND 
WRITEC 10 , 2000 , REC“6)UAVEL 
WRITECIO, 2000, REC“7)UADIR 
WRITECIO, 2000, REC=8)UATEMP 
WRITECIO,  2000,  REC-=9)UAHT 
WRITECIO, 2000, REC“10)LAPSE 
DO  40,  1-1,5 

WRITECIO, 2000 . REC-C 10+ I ) )TPRESC I ) 
WRITECIO , 2000, REC-C 15+1 ))TTEMPC I ) 
WRITECIO, 2000, REC-C20+I))TDIRC I) 
WRITEC 10, 2000, REC-C25+I))TVELC I) 


40 

CONTINUE 

2 

FORMATCAl) 

3 

FORMATCIX, ' 

INPUT  BEYOND  LIMITS',/) 

4 

FORMATCIX, ' 

Year:  ' 

,F5.0,'  GMT') 

5 

FORMATCIX, ' 

Month:  ' 

,F3.0,'  GMT') 

6 

FORMATCIX. ' 

Day:  ' 

,F3.0,'  GMT') 

7 

FORMATCIX, ' 

Time:  ' 

,14,'  GMT') 

8 

FORMATCIX. ' 

APP:  ' 

,F4.1,'  mb') 

9 

FORMATCIX,  ' 

Suface  Wind:  ' 

,F5.1,'  k.ts') 

10 

FORMATCIX,  ’ 

Wind  Speed: 

',F5.1,'  knot 

11 

FORMATCIX, ' 

Wind  Direction; 

',F5.1,'  "') 

12 

FORMATCIX, ' 

Temperature : 

',F5.1,'  “C') 

13 

FORMATCIX, ' 

Height: 

' , F7 . 1 , '  mete 

14 

FORMATCIX, ' 

Pressure ; 

',F6.1,'  mb') 

15 

FORMATCIX, ' 

Temperature ; 

',F4.1,'  “F') 

1000  F0RMATCI4) 
2000  F0RMATCF7.1) 
END 


SUBROUTINE  CLEAR 


o  o  o 


lOA 


********************************* 

THIS  S/R  CLEARS  THE  SCREEN 

********************************* 

DO  10,  1-1,30 
WRITEC*,*D'  ' 

10  CONTINUE 
RETURN 
END 


non  cj  noo  nno  o  n  o  r»  o 
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PROGRAM  WINDMAKR 

THIS  PROGRAM  READS  DIRECT  ACCESS  FILES  CREATED  BY 
'FILEMAKR'  AND  CREATES  INPUT  FILES  FOR  THE  WOCSS  WIND 
MODEL. 

A******************************** 


INTEGER  CDAY,  lYEAR,  ITIME,  JTC5D 
REAL  YEAR,  TIME,  XSC5),  YSC5),  LAPSE 
REAL  TPRESC5D  ,  TTEMPC5D  ,  TDIRCS!) .  TVELC5) 
REAL  UAHT,  UATEMP,  UADIR,  UAVEL 
OPENCUNIT  “  10, 

+  FILE  =  'MAINDATA.DAT', 

+  ACCESS  “  'DIRECT', 

+  STATUS  =  'OLD', 

+  FORM  =  'FORMATTED', 

+  RECL  =  10) 

READC10,2000,REC=13  YEAR 
READCIO, 1000,REC=2)  CDAY 
READC10,2000,REC«3:)  TIME 
READC10,2000,REC=63  UAVEL 
READC10,2000,REC*=7D  UADIR 
READC10,2000,REC-83  UATEMP 
READC10,2000,REC”9)  UAHT 
READC10,2000,REC“10)LAPSE 
DO  10,  1-1,5 

READC 10 , 2000 ,  REC-C10+  O^TPRESC  I  :> 
READC10,2000,REC-C15+n)TTEMPCO 
READC10,2000,REC=C20+n)TDIRCn 
READC  10, 2000,  REC-C25+n)TVELC  n 
CONTINUE 

OPENCll.FILE-'WINDS.DAT' , STATUS- '  NEW'  ) 


CALCULATE  DATE  AND  TIME 


IYEAR-NINTCYEAR-lOOMNTCYEAR/100)) 

IYEAR-IYEAR*1000+CDAY 

ITIME-NINTCTIME) 

WRITECll, 100)  lYEAR,  ITIME 

LIST  STATIONS 

DATA 

T ,  XS ,  YS/1 , 2 , 3  ,  <4 , 5 , 539 .6,521.5, 547 . 8 , 539 . 5 ,  557 . 1 , 3898 . 6 , 
+3895.8,3927.5,3920.2,3901.7/ 

DO  20, 1-1,5 

WRITEC11.150)  JTCI),XSCI),YSCI) 

20  CONTINUE 

ENTER  TOWER  DATA 

DO  30,  1-1,5 
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CONVERT  TTEMP  TO  CELCIUS  AND  TVEL  TO  M/S 

*ft*ft*ftft**********ft*ft***Aft**A***** 

TTEMPC I )*5 . /9 . *  CTTEMPC I D-32 . D 
TVELCO-TVELCn*  .514791 

*A*«AAAAA*AA*AAAAA*A*A*A*A**A*A*A 

WRITEC 1 1 , 200)TPRESC I ) , TTEMPC  O , TD I RC  O . TVELC I D 
30  CONTINUE 

tttt****************************** 

ASSIGN  850  WIND  TO  2000  d 

a******************************** 


UAVEL-UAVEL* .514791 

WRITECll , 250)2.1026. , TDIRC2). TVELC2) , 2000. , UADIR, UAVEL 
C  aaaaaaaa.aaa  CALCULATE  SFC  -  850  LAPSE  RATE  **»•*»*»** 

DZ-1026.-UAHT 

DTDZ-CTTEMPC2)-UATEMP)/CDZ) 

Q  aaaaaaa.aaa.aa.aaaa*  EXTRAPOLATE  TO  2000  m  *********** 

DH-2000.-UAHT 
T2-UATEMP+CDH*DTDZ) 

P2-850.-CDH* .1) 

0PENC12,FILE='TEMPS.DAT' , STATUS*' NEW' ) 

WRITEC12,300)2,2 

WRITEC 12 , 350)1026 . , TTEMPC2) ,  TPRESC2 ) 

WRITEC12. 350)2000. ,T2.P2 

11  F0RMATCA80) 

12  FORMATCAl) 

100  F0RMATCI5, IX, 14) 

150  F0RMATCI1.1X,F6.1,1X,F7.1) 

200  F0RMATCF7. 1, 1X,F5.1, 1X,F6. 1, 1X,F5. 1) 

250  F0RMATCI1,2C/.F6.0, IX, F5. 0, IX, F6 . 1)) 

300  FORMATCIl. IX, II) 

350  F0RMATCF6.0, 1X,F5. 1, 1X,F7. 1) 

1000  F0RMATCI4) 

2000  F0RMATCF7.1) 

END 


non  onwon  onn  non 
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PROGRAM  INDXMAKR 

**«*******«A*******«**il*«*******il 

PROGRAM  TO  CALCULATE  RAW  LTI 
********************************* 

REAL  MACRO . W I NDC61 ) , R0UGHC61 D , MAXGRI D 
REAL  MROUGH,  LAPSE,  INDEXC6I:) 

INTEGER  WRGRIDC61) 

OPENCIO.  FILE-' WINDVEL.DAT' . STATUS-' OLD' ) 
OPENC 11,  FI LE- ' RUFF I LE . DAT ' , STATUS- ' OLD ' D 
OPENC12,  FILE-' WR. DAT' , STATUS-' NEW' ) 
OPENC13,  FILE-' RAWINDEX.DAT' , STATUS-' NEW' ) 
OPEN  CUNIT  -  lA, 


+ 

FILE  - 

'MAINDATA.DAT' , 

+ 

ACCESS 

-  'DIRECT' , 

+ 

STATUS 

-  'OLD', 

+ 

FORM  - 

'FORMATTED' , 

+ 

RECL  - 

10) 

CALCULATE  MACRO  BTI 

MROUGH  -  5905.5/100. 

READC14,2000,REC-A)APP 

READC1A,2000,REC-5)BYSWIND 

READC1A,2000,REC-10)LAPSE 

MACRO  -  MROUGH  +  BYSWIND  +  LAPSE  +  APP 

*t1tt*t*****li*iLit*»***»t*tt*t***t*t 

READ  DATA  FROM  WIND  FILE  AND  ROUGHNESS  FILE,  ADD 
OUGHNESS 

TO  WIND  AND  FIND  MAX  VALUE 
********************************* 

MAXGRID  »  0 
DO  10  J-1,61 

READ  CIO, *)CWINDCI).  1-1,61) 

READ  Cll,  *)CROUGHCI),  1-1,61) 

DO  31  1-1,61 

WRGRIDCI)  =  NINTCWINDCI)  +  CROUGHC D/lOO)) 

IF  CWRGRIDCD.GT. MAXGRID)  MAXGRID  -  WRGRIDCI) 
31  CONTINUE 

WRITEC12, 1000)CWRGRIDCI),  1-1,61) 

10  CONTINUE 
CL0SEC12) 


CALCULATE  TURBULENCE  INDEX 


0PENC12,  FILE-'WR.DAT' ,STATUS='OLD' ) 

DO  AO  3  -  1,61 

READC12, •)CWRGRIDCI). 1=1,61) 

DO  Al,  1-1,61 

INDEXCI)  -CCWRGRIDCD/MAXGRID)  *  MACRO) 
Al  CONTINUE 

WRITEC13,A000)CINDEXCI), 1-1,61) 


L. 


AO  CONTINUE 

CL0SEC12, STATUS-' DELETE' ) 
1000  F0RMATC61C1X, lADD 
AOOO  F0RMATC61C1X,F7.1DD 
2000  F0RMATCF7.1D 
END 
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PROGRAM  WINTER? 


C 

C 

GRID. 

C 

C 


C 

C 

C 
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C 

C 

C 


********************************* 

THIS  PROGRAM  INTERPOLATES  2  KM  WINDS  TO  FILL  A  1  KM 

IT  ALSO  CALCULATES  W  COMPONENTS  FOR  THE  1  KM  GRID. 

********************************* 

REAL  VELC31,31D,DIRC31,31),UC61,61), VC61,61) 

REAL  SC61D,DC61D,DZDXC61),DZDYC61D,WC61D 

DEG2 RAD-1. /57. 295 

RAD2DEG-57.295 

OPENCIO.FILE-'SPD.DAT' ,STATUS='OLD'  ) 

OPENCll.FILE-’DIR.DAT' , STATUS-' OLD'  ) 

********************************* 

READ  2  KM  DATA 

********************************* 

DO  10,  J-31,1,-1 
READCIO, “^CVELCI , JD, 1-1,31) 

READCll. *)CDIRCI, 3), I“l,31) 

CONTINUE 

CLOSECIO) 

CLOSECll) 

********************************* 

CONVERT  TO  U  AND  V  COMPONENTS 

************lk*************t****** 

DO  20,  J-1,61,2 
DO  21,  1-1,61,2 
II-CI+l)/2 
JJ»CJ+l)/2 

UC I , J)— VELC 1 1 , JJ)*SINCDIRC 1 1 , JJ)*DEG2RAD) 

VCI . J)— VELCII , JJ)‘COSCDIRCI I, JJ)*DEG2RAD) 
CONTINUE 
CONTINUE 

********************************* 

INTERPOLATE  COMPONENTS 
********************************* 


DO  30,  J-2,60,2 
DO  31,  1-2,60,2 

UCI, J)=CUCI-1, J-1)+UCI+1,3-1)+UCI+1, J+1)+UCI- 

1, J+l))/4. 

vci,  J)-=cvci-i.  J-i)+vci +1,3-1  )+vci+i,a+i)+vci- 

l,3+l))/4. 

31  CONTINUE 
30  CONTINUE 

DO  60,  3-1,61,2 
DO  61,  1-2,60,2 

IFC3.EQ.  DTHEN 

UCI ,3)=CUCI-1,3)+UCI+1,3)+UCI ,3+l))/3. 
VCI,3)-CVCI-1,3)+VCI+1,3)+VCI ,3+l))/3. 
ELSEIFC3.EQ.61)THEN 
UCI,3)-CUCI-l,3)+UCI+l,3)+UCI,3-l))/3. 
VCI,3)-CVCI-l,3)+VCI+1.3)+VCI,3-l))/3. 
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ELSE 

UCI , JD=CUCI-1. JD+UCI+l. J)+UCI, J+ID+UCI. J-l)5/4. 

vci,  JD*cvci-i,  JD+vci+1,  j)+vci,  J+i:)+vci.  J-i)D/4. 

END  IF 

61  CONTINUE 
60  CONTINUE 

DO  70,  J=2.60,2 
DO  71,  1-1,61,2 

IFCI .EQ. IDTHEN 

UCI, J)*CUCI+1, J)+UCI, J-lD+UCI,0+lD)/3. 

VCI, JD«CVCI+1, J)+VCI, J-1)+VCI, J+l))/3. 

ELSEIFCI .EQ.61DTHEN 
UCI, JD-CUCI-l.JD+UCI, J-1)+UCI, J+l)D/3. 

VCI, JD-CVCI-1, J)+VCI, J-1)+VCI, J+l))/3. 

ELSE 

uci,  JD»cuci-i,  J)+uci+i.  j)+uci,  J+iD+uci ,  j-id:)/4. 
vci , JD*CVCI-1, J)+VCI+1, J)+VCI, J+ID+VCI, J-l))/4. 
END  IF 

71  CONTINUE 
70  CONTINUE 

TRANSLATE  BACK  TO  METEOROLOGICAL  WINDS  AND  WRITE  TO 

LES 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

0PENC15, FILE-' WINDVEL.DAT' , STATUS-' NEW'  ) 

OPENC 1 6 ,  F I  LE= '  W I NDD I R .  DAT '  ,  STATUS- '  NEW '  ]) 

DO  50,  J-61,1,-1 
DO  51,  1-1,61 

IFCUCI ,  J:).EQ.0.0.  AND.  VCI ,  J).EQ.O.O)THEN 
SC  0=0.0 
DCO=0.0 
ELSE 

SCO=SQRTCUCI,  JD**2+VCI ,  J)**2) 

O-AM0DC540.+ATAN2CUCI ,  J),  VCI ,  J:))*RAD2DEG,360.  ) 

SCO=SCO*  .  194254 
END  IF 

51  CONTINUE 

WRITEC15,1000DCSCO.  1  =  1,60 
WRITEC16,  1000:)CDCO.  1  =  1,60 
50  CONTINUE 
CL0SEC15D 
CL0SEC16) 

**********************tk********** 

CALCULATE  W'S  AND  WRITE  TO  FILE 
*••******«***•*•*•«************** 

0PENC12,FILE='X.DAT' .STATUS-' OLD' ) 

0PENC13, FILE-' Y. DAT' , STATUS-' OLD'  ) 
0PENC14,FILE='W.DAT' . STATUS-' NEW'  D 
DO  40,  J-61,1,-1 
READC12,  OCDZDXCO,  1  =  1,60 
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READC13,  *)CDZDYCn.  I“1,61D 

A******************************** 

TAKE  DOT  PRODUCT  OF  WIND  AND  TERRAIN,  THEN  CONVERT  TO 

KTS 

Q  AA***«**AA**A*AAAAAAA*AAAA*AAAAAA 

DO  41,  1=1,61 

WCI)*CUCI, J3‘DZDXCn+VCI,JD *DZDYC I .000194254 
41  CONTINUE 

WRITEC14,  lOOO^CWCn,  1  =  1, 6i:) 


40  CONTINUE 
1000  F0RMATC61C1X,F7.1DD 
END 


non  oon  ooo  ooon 


112 


PROGRAM  INDXMOD 

PROGRAM  TO  MODIFY  LTI  BASED  ON  WAKE  AND  PI  REPS 

OUTPUT  IS  FOR  CAT  1  AIRCRAFT 

********************************* 

REAL  YEAR,  MYEAR,  MTIME,  TIME,  YE,  TI,  INDEXC61,61) 
REAL  YClO),  XCIO),  TCIOD,  WC61).  XX,  YY,  FACTOR 
INTEGER  CDAY,CA,MCDAY, CHECK, ACCAT, INC lOD , ACC  10) 
INTEGER  PIREP, TERRTYPEC61, 61) , AT, TME 

CHARACTER  INTENS,  ID 
********************************* 

LOAD  INDEX  ARRAY 

********************************* 

OPENC 13 , F I LE= ' RAW I NDEX . DAT ' , STATUS= ' OLD ' ) 

DO  AO,  J=61, 1,-1 
READC13, *)CINDEXCI , J), 1=1,61) 

CONTINUE 

********************************* 

ADD  20  FOR  NEGATIVE  W 
********************************* 

0PENC1A,FILE=' W.DAT' , STATUS=' OLD'  ) 

MARKER=0 
DO  50,  J=61,l,-1 
READC14, AOOO)CWCI), 1=1,61) 

DO  51,  1=1,61 

IFCWCI).LT.-3.)THEN 
INDEXCI. J)=INDEXCI. J)+20. 

MARKER=MARKER+1 
END  IF 
51  CONTINUE 
50  CONTINUE 

IFCMARKER.EQ.O)THEN 
WRITEC * , * ) '  No  wake  modifications.' 

END  IF 

**********'»********************** 

FIGURE  SEARCH  TIMES  FOR  PI REPS 
********************************* 

OPEN  CUNIT  =  10, 

+  FILE  =  ’MAINDATA.DAT', 

+  ACCESS  =  'DIRECT', 

+  STATUS  =  'OLD', 

+  FORM  =  'FORMATTED', 

+  RECL  =  10) 

READC10,2000,REC=1)YEAR 
READCIO, 1000,REC=2)CDAY 
READCIO, 1000, REC=3)TME 
TIME=FLOATCTME) 

MTIME=TME-100 
MCDAY=CDAY 
MYEAR=YEAR 

IF  CTME  .LT.  100)  THEN 
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MTIME=2300+TME 
MCDAY=CDAY-1 
IF  CCDAY.EQ.ID  THEN 

YRES  =  AMODCCYEAR+1 0,4.0) 

GRES  =  AMODCCYEAR+1. ). 100.0) 

MCDAY=365 

I F  C YRES . EQ . 0 . 0 . AND . GRES . NE . 0 . 0)MGDAY=366 
MYEAR=YEAR-1. 

END  IF 
END  IF 

LOOK  FOR  PI REPS 

0PENCUNIT=11, 

+  FILE='PIREP.DAT' , 

+  AGGESS=' DIRECT' , 

+  STATUS='OLD' , 

+  FORM=' FORMATTED' . 

+  RECL=50) 

ICT=0 

DO  10,  1=1,999 
IND=I-1 

READCll, 101,REC=I ,END=11)YE,CA,TI ,XX,YY, ID, AT 
IFCYE.EQ.9999. )G0  TO  10 

IFCCYE.LT.MYEAR).OR. CCA . LT . MCDAY) . OR. CTI . LT . MTIME))THEN 
ICT=1+ICT 
GO  TO  10 
END  IF 

IND=IND-ICT 

IFCCID.EQ. 'L' ).0R. CID.EQ. ' 1 ' )) INC IND)=1 
IFCCID.EQ. 'M').OR.CID.EQ. 'm'))INCIND)=2 
IFCCID.EQ. 'S' ). OR. CID.EQ. 's' ))INCIND)=3 
IFCCID.EQ. 'X' ). OR. CID.EQ. 'x' ))INCIND)=4 
INCIND)=INCIND)+CAT-1) 

TCIND)=TI 

XCIND)=XX 

YCIND)=YY 

10  CONTINUE 

11  IFCIND.EQ.O)  THEN 
WRITEC* , *)'There  are  no  current  PIREPS.' 
0PENC15,FILE=' NEWINDEX.DAT' . STATUS=' NEW' ) 

DO  70,  3=61, 1,-1 

WRITEC15,4000)CINDEXCI, 3), 1=1,61) 

CONTINUE 
GO  TO  9999 
ELSE 

IND=IND-ICT-1 

IFCCID.EQ. 'L' ). OR. CID.EQ. ' 1 ' )) INC IND)=1 
IFCCID. EQ. 'M' ).0R. CID.EQ. 'm' ))INCIND)=2 
IFCCID.EQ. 'S' ). OR. CID.EQ. ' s ' )) INC IND)=3 


70 
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IFCCID.EQ. 'X' ).OR.CID.EQ. 'x' )DINCIND)=4 
INCIND)=INCIND>+CAT-1D 
TCINDD=TI 
XCIND)=XX 
YCIND)-YY 
END  IF 

********************************* 

READ  TERRTYPE  FILE 
********************************* 

0PENC12,  FILE*' TERRTYPE. DAT' , STATUS*' OLD'  D 
DO  18,  1=1,61 

READC12,*D  CTERRTYPECI, JD,  J=l,61) 

18  CONTINUE 

********************************* 

MODIFY  LTI 

********************************* 

DO  20,  N=l, IND 

********************************* 

FIND  NEAREST  GRIDPOINT 
********************************* 

X 

********************************* 

XMIN=506. 1 

XX=XCN)-XMIN 

IX=NINTCXXD 

********************************* 

Y 

A******************************** 

YMIN=3887.0 

YY=YCND-YMIN 

IY=NINTCYY) 

CHECK  FOR  PI REP  >  INDEX 

PIREP  ASSIGNED  MAXIMUM  BTI  VALUE  FOR  CATEGORY 

*************A*AA***A*****»**AA** 

PIREP=59. 

DO  30,  L=l, INCND 

PIREP=FL0ATCL*10D+PIREP 
30  CONTINUE 

IFCINDEXCIX, IY3.GE.PIREP)G0  TO  20 
FACTOR=PIREP/INDEXCIX, lYD 

INDEXCIX, IY)=PIREP 
********************************* 

MODIFY  INDEX  AT  LIKE  TERRAIN 
********************************* 

DO  21,  1=1,61 

DO  22,  J=l,61 

IFCI .EQ. IX.AND. J.EQ. lYDGO  TO  22 
IFCTERRTYPECI ,  J:).EQ.TERRTYPECIX,  IY^DTHEN 
INDEXCI , J)=CINDEXCI . J)*FACTOR) 

END  IF 


115 


22  CONTINUE 

21  CONTINUE 
20  CONTINUE 

0PENC15,FILE=' NEWINDEX.DAT' , STATUS=' NEW' D 
DO  60.  J=61,l.-1 

WRITEC15,4000)CINDEXCI. JD. 1=1, 61D 
60  CONTINUE 

101  F0RMATC1X,F6.0.2X, 13 , F6 . 0, 2X, 2CF11 . 1 , 2X) , A1 , 2X , ID 

1000  F0RMATCIA3 

2000  F0RMATCF7.1D 

3000  F0RMATC61CIADD 

4000  F0RMATC61C1X,F7.  DD 

9999  END 


no  o  on  o  on  on  «» o  o  o  o  o  o  o 
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PROGRAM  LT I GRAPH 

MODIFIED  TO  28x28  GRID  AROUND  19  JULY  1988 
MODIFIED  28  JULY  1988 

SO  THAT  DATA  FILES  RESIDE  IN  DATA  SUBDIRECTORY 

MODIFIED  FOR  61x61  GRID  3  FEB  1989  BY  N.  SMITH 


PROGRAM  TO  PLOT  TURBULENCE  INDEX 


NOTRUNCATE 

dimension  zC61 , 61 J , vertC61 , 61 J 
REAL  ZTERMI , 2TERMA,  V_MIN,  V_MAX 

INTEGER  ITGRMI.  JTGRMI ,  ITGRMA,  JTGRMA,  IVGRMI,  JVGRMI 
INTEGER  IVGRMA,  JVGRMA,  LEVEL 
COMMON  ZTERM I , I TGRM I . JTGRM I . 

1  ZTERMA , I TGRMA , JTGRMA , 

2  V_MIN. IVGRMI, JVGRMI. 

3  V_M AX . I VGRM A , JVGRMA , 

4  TIGRMI .TJGRMI.TIGRMA.TJGRMA, 

5  VIGRMI , VJGRMI, VIGRMA, VJGRMA 


-  DO  THE  TERRAIN  CONTOURS  - 

openClO, f 1 le=' NTCIOOO. PLT' ) 

READING  THE  TERRAIN  FILE 
wrlteC6, Reading  the  terrain  file...' 
readClO,*)  swx, swy, numx, numy , d 
do  100  J=61,l,-1 

readC10,*J  CzCI , JJ , 1=1 , 61) 

100  continue 

Z  IS  THE  FIRST  FILE  TO  BE  WRITTEN  TO  THE  DATA  FILE 
f I le*l 

call  contourCz , f 1 le) 


-  DO  THE  VERTICAL  WIND  FIELD  CONTOURS  - 

wr 1 teC6, 90) 

90  formate// , IX, ' Reading  the  turbulence  file...') 
0PENC22,FILE«' NEWINDEX.DAT' , STATUS-' OLD'  ) 

DO  1010  J=61,l,-1 
READC22,4000)  CVERTC I , J) , 1=1 . 61) 

4000  F0RMATC61C1X,F7.1)) 

1010  CONTINUE 

VERT  IS  THE  SECOND  FILE  TO  BE  WRITTEN  TO  THE  DATA  FILE 
f lle=2 

call  contourCVert , f l le) 

CALL  LABEL 
end 


SUBROUTINE  CONTOURCZ, fl le) 


n  ac  n 
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REAL  ZC61,61) 

DIMENSION  XXC5:),YYC5D 
DIMENSION  XC6i:),YC61) 

CHARACTER* 15  FILENAME 
CHARACTER* 10  ACCES 
CHARACTER*!  AELSTP 

COMMON  ZTERMI. ITGRMI, JTGRMI, 

1  ZTERMA , I TGRMA , JTGRM A , 

2  V_MIN. IVGRMI, JVGRMI , 

3  V_MAX . I VGRMA . JVGRMA , 

4  TIGRMI,TJGRMI,TIGRMA,TJGRMA, 

5  VIGRMI.VJGRMI.VIGRMA.VJGRMA 

X9-1.0 
Y9*1.0 
DZ— .021 
SCL=30./60. 

FILENAME"' PLOT. DAT' 

IFCFILE  .EQ.  1.)  THEN 

C - SET  UP  FILE  FOR  TERRAIN  DATA - 

ACCES  "  'SEQUENTIAL' 

NSTART"! 

NP0INTS"61 
ELSE 

- SET  UP  FILE  FOR  THE  SECOND  SET  OF  DATA,  THE  VERTICAL 

IND  DATA- 


POSITION  FILE  POINTER  TO  END  OF  FILE - 

ACCES  =  'SEQUENTIAL' 

10  READCLUOUT, 'CIXD' ,END=20) 

GOTO  10 

20  NSTART"! 

NP0INTS"61 

C - BACKSPACE  TO  WRITE  OVER  THE  END-OF-FILE  MARKER 

BACKSPACE  8 
END  IF 

C  OPEN  PLOT  FILE  TO  BE  WRITTEN 


LUOUT-8 

OPENCUNIT=LUOUT,FILE=FILENAME, ACCESS=ACCES) 
IFCFILE  .EQ.  1.)  WRITECLU0UT,45)' AF;  IN;  ' 
NP0INTS"61 

DO  90  I"NSTART,NPOINTS 
XCn"FLOATCI-lD 
YCO=FLOATCI-1D 
90  CONTINUE 
140  F0RMATC1X,61F5.0,/) 

C  INITIALIZE  PLOTTER 

WRITECLU0UT,85D  'IN;' 

WRITECLU0UT,55)  ' IPIOOO, 1000,7000,7000;  ' 
WRITECLU0UT,65)  ’ SCO, 30, 0, 30; ' 

45  F0RMATCA6) 

55  F0RMATCA22) 
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65  F0RMATCA12D 
75  FORMATCAAD 
85  FORMATCA3) 

95  F0RMATCA7) 

IFCflle  .eq.  ID  CALL  FRAMECLUOUTD 

C - 

WRITECLU0UT,85D'SP6; ' 

C  If  not  the  first  file,  choose  pen  2 

IFCflle  .ne.  lO  WRITECLUOUT,  75D  'SP2;' 
ZMTN-9999999 
2MAX»-99999 
IFCFILE  .EQ.  ID  THEN 
DO  150  I»l,61 
DO  150  J*l,61 
2MIN=AMIN1CZMIN,ZCI, JDD 
ZMAX-AMAXICZMAX. ZCI.JDD 
150  CONTINUE 


DO  160  1-1,61 
DO  160  J=l,61 
IFCZCI.JD  .EQ.  ZMIND  THEN 
IMIN=I 
JMIN-J 
END  IF 

IFCZCI,JD  .EQ.  ZMAXD  THEN 
IMAX-I 
JMAX-J 
END  IF 

160  CONTINUE 

TIGRMI  -  FLOATCIMIN-lD+506. 1 
TJGRMI  -  FLOATCJMIN-lD+3887.0 
WRITE  C6,280DZMIN, TIGRMI, TJGRMI 
ZTERMI  -  ZMIN 
ITGRMI  -  IMIN-1 
JTGRMI  -  JMIN-1 
TIGRMA  -  FLOATCIMAX-lD+506. 1 
TJGRMA  -  FLOATCJMAX-lD+3887.0 
WRITE  C6,285DZMAX, TIGRMA, TJGRMA 
ZTERMA  -  ZMAX 
ITGRMA  -  IMAX-1 
JTGRMA  -  JMAX-1 
END  IF 

IFCFILE  .EQ.  2D  THEN 
DO  151  1=1,61 
DO  151  J=l,61 

ZMIN=AMIN1CZMIN,ZCI. JDD 
ZMAX-AMAX1CZMAX,ZCI , JDD 
151  CONTINUE 

DO  161  1=1,61 
DO  161  J-1,61 


o  n 
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IFCZCI.JD  .EQ.  ZMIN3  THEN 
IMIN*I 
JMIN-J 
END  IF 

IFCZCI.JD  .EQ.  ZMAXD  THEN 
IMAX-I 
JMAX*J 
END  IF 

161  CONTINUE 

VIGRMI  =  FLOATCIMIN-lD+506. 1 

VJGRMI  »  FL0ATCJMIN-l)+3887.0 

WRITE  C6,290:)ZMIN,  VIGRMI,  VJGRMI 

V_MIN  -  ZMIN 

IVGRMI  *  IMIN-1 

JVGRMI  »  JMIN-1 

VIGRMA  =  FLOATCIMAX-lD+506. 1 

VJGRMA  -  FL0ATCaMAX-l)+3887.0 

WRITE  C6,295DZMAX, VIGRMA, VJGRMA 

V_MAX  =  ZMAX 

IVGRMA  =  IMAX-1 

JVGRMA  =  JMAX-1 

END  IF 

280  FORMATC/,'  Min  ',flO.A,’  m  at  C ' , F5 . 1 . ' , ' , F6 . 1 , ' J ) 
285  FORMATC'  Max  ',fl0.4,'  m  at  C ' , F5 . 1 , ' , ' , F6 . 1 , ' ) ' ) 

290  FORMATC/,'  Min  ',fl0.4,'  at  C ' , F5 . 1 , ' , ' , F6 . 1 ,  '  ) '  J 
295  FORMATC'  Max  ’,fl0.4,'  at  C ' . F5 . 1 , ' , ' , F6 . 1 , ' ) ' J 

Default  value  for  elstp  to  generate  10  contour  lines 
ELSTP-CZMAX-ZMIN)/10. 

WRITEC6, 175J  ELSTP 
175  FORMATCAl) 

IFCFILE  .EQ.  IJ  THEN 


ELSTP-200 

c 

WRITEC6, * J 

9 

You 

may  choose: 

9 

c 

WRITEC6, 

9 

a 

100 

meters ' 

c 

WRITEC6. * J 

9 

b 

200 

meters' 

c 

WRITEC6, • J 

9 

c 

500 

meters ' 

c 

WRITEC6, * J 

9 

d 

1000 

meters' 

c 

WRITEC6, • J 

9 

e 

no  contours' 

c 

c 

WRITEC6, * J 
WRITEC6, * J 

9 

for 

contour  Interval. 

9 

c 

c 

WRITEC6, * J 
EEADC5, 175J 

'  Input 
AELSTP 

a , b, c , d,  or 

e; 

c 

IFCAELSTP  . 

EQ. 

'  A' 

.OR.  AELSTP 

.EQ. 

'a'J  ELSTP=100. 

c 

IFCAELSTP  . 

EQ. 

'B' 

.OR.  AELSTP 

.EQ. 

'b'J  ELSTP=200. 

c 

IFCAELSTP  . 

EQ. 

'C' 

.OR.  AELSTP 

.EQ. 

'c'J  ELSTP=500. 

c 

IFCAELSTP  . 

EQ. 

'  D' 

.OR.  AELSTP 

.EQ. 

'd'J  ELSTP=1000. 

c 

IFCAELSTP  . 
ELSE 

EQ. 

'E' 

.OR.  AELSTP 

.EQ. 

'e'J  ELSTP=900000 

ELSTP=10. 

IFCMAXCABSCZMINJ,ABSCZMAX)J  .GE.120.0J  THEN 
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WRITEC6,  *  ]) '  -»»»»»«=»«■*•*•===-»==»===>==========*========*=■!==' 

WRITEC6, *D 
WRITEC6, *) 

WRITEC6, *D 
WRITEC6, 

WRITEC6. 178D 

178  FORMATC'  - CAUTION:  SOME  EXTREME  TURBULENCE - 

WRITEC6, 

WRITEC6, 

WRITEC6, 

WRITEC6, 

^WRITECe, 

PAUSE'  Press  Enter  to  continue...’ 

END  IF 

IFCZMAX-ZMIN  . LT .  1.)  THEN 
WRITEC6,  *:) 

WRITEC6, 

WRITEC6, *) 

WRITEC6, 

WRITEC6. O 

WRITEC6,  *)'  *========*=======»=-=====«======«*===*=*■«======' 

WRITEC6, 

WRITEC6, 177) 

177  FORMATC'  Range  of  turbulence  index  Is  less 

than'  , 

&  '  1  .  '  ) 

WRITEC6,*)C’  The  turbulence  index  Is  uniform.') 
WRITEC6, •) 

WRITEC6,  *  )'  »»**=«==»===»==»=»=«==»==>====*==============' 

WRITEC6, ‘) 

WRITEC6, *) 

WRITEC6, •) 

WRITEC6, *) 

PAUSE'  Press  Enter  to  continue...' 

CLOSECUNIT-LUOUT, STATUS*' DELETE' ) 

GOTO  1020 
END  IF 

WRITEC6, 176)ELSTP 

176  FORMATC’  Contour  level  for  turbulence  Index  Is  ',F5.2, 
&  '  .  '  ) 

ENDIF 

c  IFCELSTP  .LT.  CZMAX-ZMIN)/30. )  ELSTP=CZMAX-ZMIN)/10 . 
MI-AINTCZMIN/ELSTP) 

MA-AINTCZMAX/ELSTP) 

WRITEC6,*)  'Generating  Plot  file  from  data  files...  ' 


121 


DO  1010  K-MI,MA 

IFCFILE.EQ.2)  THEN 

IFCK.EQ.6.0R.K.EQ.7.0R.K.EQ.9.0R.K.EQ. lO.OR.K.EQ. 12)THEN 
GO  TO  A44 
ELSE 

GO  TO  1010 
END  IF 

444  ENDIF 

F-FLOATCK)‘ELSTP+DZ 
IFCK.eq.lOD  THEN 

IFCf  lle.EQ.2:)WRITECLU0UT,75D'SP4;  ' 

ENDIF 

IFCK.eq.7D  THEN 

IFCf 1 le . EQ. 2DWRITECLU0UT, 755' SP3 ; ' 

ENDIF 

IFCK.ge.l2D  THEN 

IFCf 1 le . EQ. 2)WRITECLU0UT, 75) ' SP5; ' 

ENDIF 

IFCK  .le.  MI .or. k. It .7)  THEN 
IFCflle  .EQ.  2)WRITECLU0UT,75)'SP2; ' 

IFCf lie  .EQ.  1)WRITECLU0UT,75)'SP6; ' 

ENDIF 

IFCK.eq.9)  THEN 

IFCflle  .EQ.  2)WRITECLU0UT,75)'SP8j  ' 

ENDIF 

WRITECLU0UT,85)  'PU;' 

DO  1000  1-1,60 
DO  1000  J-1,60 
FO-ZCI, J)-F 
Fl-ZCI+1. J)-F 
F2-ZCI, J+l)-F 
F3-ZCI+1, J+ID-F 
IC-1 
G-F0*F1 

IFCG.LE.0.0)  THEN 

XXCIC)-CCXCn*Fl-CXCI)+X9)*F0)/CFl-FQ))*SCL 
YYCIC)-CYCJ))*SCL 
IC-IC+1 
END  IF 
G-F2*F0 

IFCG.LE.0.0)  THEN 
XXCIC)-CXCI))*SCL 

YYCIC)-CCYCJ)*F2-CYCJ)+Y9)*F0)/CF2-F0))*SCL 
IC-IC+1 
END  IF 
G-F3*F1 

IFCG.LE.0.0)  THEN 

XXC IC)-CXC I )+X9) ‘SCL 

YYCIC)-CCYCJ)*F3-CYCJ)+Y9)*F1)/CF3-F1))*SCL 

IC-IC+1 


L. 
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END  IF 
G“F3*F2 

IFCG.LE.O.OD  THEN 

XXCICD-CCXCn*F3-CXCi:)+X9:)*F2D/CF3-F2D:)*SCL 
YYCIC)-CYCJ:)+Y9)*SCL 
IC-IC+1 
END  IF 

IFCIC.EQ.ID  GO  TO  1000 
IFCIC.EQ.5D  GO  TO  2000 
CALL  F0RMTCXXC1D,YYC1D,LU0UT) 

WRITECLU0UT,85)  'PD;' 

CALL  F0RMTCXXC2:),YYC2D.LU0UT) 

WRITECLU0UT.85:)  'PU;' 

GO  TO  1000 

2000  IFCFO.LE.O.  AND.F3.LE.0.0:)  THEN 
DO  2001  M«l,2 

N-H+2 

CALL  FORMTCXXCMD,YYCM).LUOUT) 

WRITECLU0UT,85D  'PD;' 

CALL  FORMTCXXCND.YYCN),LUOUTD 
WRITECLU0UT,85D  'PU;' 

2001  CONTINUE 
ELSE 

DO  2002  M«l,3,2 
N-M+1 

CALL  FORMTCXXCM),YYCM).LUOUT) 

WRITECLU0UT,85D  'PD;' 

CALL  FORMTCXXCND.YYCND.LUOUT) 

WRITECLU0UT,85D  'PU;' 

2002  CONTINUE 
END  IF 

1000  CONTINUE 

IFCK  .EQ.  0  .OR.  K  .LE.  MI  .OR.  K  .EQ.  MA)  THEN 
IFCflle  .EQ.  1)  WRITECLU0UT,75)’SP6; ' 

IFCflle  .EQ.  2)  WRITECLUOUT, 75) ' SP2 ; ' 

END  IF 

1010  CONTINUE 

IFCFILE  .EQ.  2)  WRITECLUOUT, 85)' AF; ' 

1020  CONTINUE 

RETURN 
END 

C - 

SUBROUTINE  FRAMECLUOUT) 

C  DRAW  THE  THE  FRAME 

WRITECLU0UT.5) 

' SPl ; PU ; PAO, 0; PD; PAO, 30; PA30, 30; PA30, 0; PAO, 0; PU;  ' 

5  F0RMATCAA7) 

DO  100  I»l,31 

IFCI  .LT.  11)  WRITECLUOUT, 10)1-1, '0' , I-l, '0.5' 

100  IFCI  .GE.  11)  WRITECLUOUT, 20)1-1, 'O’ , I-l, '0.5' 


no  X  X  X  X  no 
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10  FORMATC'PA' . II. ' , ' , Al. ' ; ' , 

'PD' , ' } ' . 'PA' , II. ' , ' .A3. ' ; ' . 'PU' , ' ; ' ) 

20  FORMATC’ PA' . 12. ’ . ' . Al. ' ; ' . 

'PD' . ' j ' . 'PA' . 12. ' . ' . A3. ' ; ' . 'PU' . ' ; ' ) 

DO  200  1-1.31 

IFCI  .LT.  11)  WRITECLU0UT.30)'30' . I-l 
200  IFCI  .GE.  11)  WRITECLU0UT.'40)'30 

30  FORMATC'PA' .A2. '.'. II. . 

'PD' . ' ; ' . 'PA' .A4. ' , ' . II. ' ; ' . 'PU' , ' ; ') 

40  FORMATC'PA' .A2. '.'. 12. . 

'PD' . ' ; ' . 'PA' . A4. ' . ' , 12. ' ; ' . 'PU' . ' ; ' ) 

DO  300  1-31. l.-l 

IFCI  .LT.  11)  WRITECLUOUT. 50)1-1. '30' 
300  IFCI  .GE.  11)  WRITECLUOUT. 60)1-1 

50  FORMATC'PA' . II. '.' .A2. . 

'PD' . ' ; ' . 'PA' . II. ' . ' .A4. ' ; ' . 'PU' . ' ; ' ) 

60  FORMATC'PA' . 12. '.' .A2. . 

'PD' . ' ; ' . 'PA' . 12. ' . ' .A4. ' ; ' . 'PU' . ' ; ' ) 

DO  400  1=31. l.-l 

IFCI  .LT.  11)  WRITECLUOUT. 70)'0' . I-l, 
400  IFCI  .GE.  11)  WRITECLUOUT. 80)' 0' 

70  FORMATC'PA' .Al. '.'. II. . 

'  PD'  .  '  j  '  .  '  PA'  .  A3 .  '  .  '  .  1 1 .  '  ;  '  .  '  PU'  .  '  ;  '  ) 

80  FORMATC'PA' .Al. '.'. 12, . 

'PD' . ' ; ' , 'PA' ,A3, ' , ' , 12, ' ; ' , 'PU' , ' ; ' ) 

C - DONE  DRAWING  FRAME - 

CALL  MAPCLUOUT) 

RETURN 

END 


,  '29.5' , I-l 
'  ,  I-l, '29.5'  ,  I-l 


.  I-l, '29.5' 

,  '30'  , I-l. '29.5' 


'0.5' , I-l 
,  I-l, '0.5' , I-l 


SUBROUTINE  FORMTCX , Y , LUOUT) 


IF 

cx 

.LT. 

9.995 

.AND. 

Y 

.LT. 

9.995) 

WRITECLUOUT, 1100) 

IF 

cx 

.LT. 

9.995 

.AND. 

Y 

.GE. 

9.995) 

WRITECLUOUT, 1101) 

IF 

cx 

.GE. 

9.995 

.AND. 

Y 

.LT. 

9.995) 

WRITECLUOUT, 1110) 

IF 

cx 

.GE. 

9.995 

.AND. 

Y 

.GE. 

9.995) 

WRITECLUOUT,  1111) 

,Y 

1100  FORMATC'PA' , F4 . 2 , ' , ' , F4 . 2 ,  '  ;  '  ) 

1101  FORMATC'PA' , F4. 2, ' , ' , F5. 2, ' ; ' ) 

1110  FORMATC'PA' , F5 . 2 , ' , ' . F4 . 2 . ' ; ' ) 

1111  FORMATC'PA' , F5 . 2 , ' , ' , F5 . 2 ,  '  ;  '  ) 
RETURN 

END 


SUBROUTINE  LABEL 
COMMON  ZTERMI, ITGRMl, JTGRMI, 

1  ZTERMA, ITGRMA. JTGRMA, 
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2  V_MIN, IVGRMI, JVGRMI, 

3  V_MAX , I VGRMA , JVGRMA , 

4  TIGRMI.TJGRMI.TIGRMA.TJGRMA, 

5  VIGRMI,VJGRMI.VIGRMA,VJGRMA 
CHARACTER‘80  STRING 

CHARACTER* 15  FILENAME 
CHARACTER* 10  LAB 
FILENAME-' PLOT. DAT' 

LUOUT-8 

C  DRAW  BOX 

CALL  BOXCLUOUT, 1,31. , 0 . , 45 . , 0 . , 45 . , 30 . 1 , 31 . , 30 . ID 
CALL  BOXCLUOUT, 1,31. , 28 . , 45 . , 28 . , 45 . , 30 . 1 , 3 1 . , 30 . ID 
C  WRITECLUOUT , *  D ' SP8 ; PU ; PA3 1 , 0 ; PD ; PA47 , 0 ; PA47 , 30 ; PA3 1,30; 

PA31,0; ' 

C - ASL  LABEL - 

WRITECLUOUT, *D  'SP5;' 

WRITECSTRING,99D 

CALL  LBCOMMANDCLUOUT, 31 . 5,29.3, STRINGD 

C - FIRST  LABEL - 

WRITECLUOUT, *D  'SP4;' 

WRITECSTRING, lOOD 

CALL  LBCOMMANDCLUOUT, 31. 5, 27. , STRINGD 

WRITECSTRING, lOlD  ZTERMA, '  m  ' 

CALL  LBCOMMANDCLUOUT, 31. 5, 26. , STRINGD 

WRITECSTRING, 102D  TIGRMA, TJGRMA 
CALL  LBCOMMANDCLUOUT, 31. 5, 25. , STRINGD 

C - SECOND  LABEL - 

WRITECSTRING, 103D 

CALL  LBCOMMANDCLUOUT, 31 .5,23., STRINGD 

WRITECSTRING, lOlD  ZTERMI,'  m  ' 

CALL  LBCOMMANDCLUOUT, 31. 5, 22. , STRINGD 

WRITECSTRING, 102D  TIGRMI , TJGRMI 
CALL  LBCOMMANDCLUOUT, 31. 5, 21. , STRINGD 

C - THIRD  LABEL - 

WRTTFf^TRTNR  104'; 

CALL  LBCOMMANDCLUOUT, 31. 5, 19. , STRINGD 
WRITECSTRING, 105D  V_MAX, '  ' 

CALL  LBCOMMANDCLUOUT, 3 1.5, 18. , STRINGD 

WRITECSTRING, 102D  VIGRMA, VDGRMA 
CALL  LBCOMMANDCLUOUT, 31. 5, 17. , STRINGD 

C - FOURTH  LABEL - 
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WEITECSTRING. 106D 

CALL  LBC0MMANDCLU0UT,31.5. 15. , STRING) 

WRITECSTRING, 105)  V_MIN, '  ' 

CALL  LBC0MMANDCLU0UT,31.5,14. , STRING) 

WRITECSTRING, 102)  VIGRMI , VJGRMI 
CALL  LBC0MMANDCLU0UT,31.5.13. .STRING) 

C - FIFTH-NINTH  LABELS - 

WRITECLUOUT, *)  'SP6;' 

WRITECSTRING, 107) 

CALL  LBC0MMANDCLU0UT,31.5, 10. 5. STRING) 

WRITECLUOUT.*)  'SP5;' 

WRITECSTRING, 108) 

CALL  LBC0MMANDCLU0UT,31.5,9. .STRING) 

WRITECLUOUT,*)  'SP4;' 

WRITECSTRING, 109) 

CALL  LBCOMMANDCLUOUT .31.5,7.5, STRING) 

WRITECLUOUT,*)  ' SP8 ; ' 

WRITECSTRING, 110) 

CALL  LBCOMMANDCLUOUT, 31. 5. 6. .STRING) 


WRITECLUOUT,*)  ' SP3 ; ' 

WRITECSTRING,  111) 

CALL  LBCOMMANDCLUOUT, 31 .5,4 . 5, STRING) 

WRITECLUOUT,*)  'SP2;' 

WRITECSTRING, 112) 

CALL  LBCOMMANDCLUOUT, 31. 5, 3. .STRING) 

C  REWINDC22) 

C  READC22,*)  LEVEL 

C  WRITECLUOUT,*)  ' SP3 ; ' 

C  WRITECSTRING, 113)  '  LEVEL', LEVEL 

C  CALL  LBCOMMANDCLUOUT, 31. 5, 1.3, STRING) 

99  FORMATC'**  TURBULENCE  INDEX  **') 

100  FORMATC ' Terra  1 n  Maximum  of') 

101  F0RMATCF7. 1, A3) 

102  FORMATC 'at  C ' , F5 . 1 , ' , ' , F7 . 1 , ' ) ' ) 

103  FORMATC 'Terra in  Minimum  of') 

104  FORMATC ' Turb .  Index  maximum  of') 

105  F0RMATCF7.2. A5) 

106  FORMATC ' Turb .  Index  minimum  of') 

107  FORMATC ' Terra  In  Contours') 

108  FORMATC 'Extreme  Turb.  >  120') 


L 


no  on  o  on  on  o  o  no 
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109 

FORMATC ' Severe  Turb. 

100-119'  ) 

110 

FORMATC'Mod.-Sev.  Turb 

.  90-99') 

111 

FORMATC 'Mod.  Turb. 

70-89'  ) 

112 

FORMATC' Light  Turb. 

60-69'  ) 

113 

F0RMATCA17.  13) 

WRITECLUOUT, 'SPOjAFj' 

RETURN 

END 


SUBROUT I NE  LBCOMMANDCLUOUT , X , Y , STRING) 


INTEGER  LUOUT 
REAL*4  X.Y 

CHARACTER*80  STRING 


CHARACTER*!  TERMINATOR 


TERM  I NAT0R=CHARC3 ) 

WRITECLUOUT, * j  ' PU ; PA' , X, ' , ' , Y . ' ; ' , 

2  'LB' ,STRINGC1:22), TERMINATOR,  '  ;PU;  ' 

RETURN 

END 


SUBROUTINE  BOXCLUOUT, PEN. XI , Y1 , X2 , Y2 , X3 , Y3 , X4 , Y4) 

DRAW  BOX 


INTEGER  LUOUT. PEN 

REAL*4  X1.Y1,X2,Y2,X3, Y3,X4, Y4 


WRITECLUOUT. *) 
2 
3 

RETURN 

END 


'SP' .PEN, ' ;PU;PA' ,X1, ' , ' ,Y1,  '  ;PD;  '  , 
'PA' ,X2. ' , ' .Y2. ' , ' ,X3, ' , ' .Y3, ' . ' , 
X4, ' , ' ,Y4, ' , ' ,X1, ' , ' , Yl. ' ;PU; ' 


SUBROUTINE  MAPCLUOUT) 
DRAWS  MAP  OF  NTC 


INTEGER  LUOUT 
WRITECLUOUT,*)  'SP1;PU;' 
WRITECLUOUT, *)' PA6.05, .7;PD; ' 
WRITECLUOUT, * ) ' PA6 . 05, 2 . 35; ' 
WRITECLUOUT, * ) ' PA5 . 05, 2 . 35; ' 
WRITECLUOUT, * ) ' PA5 . 05, 8 . 9; ' 
WRITECLUOUT, * ) ' PA3 . 6, 8 . 9;  ' 
WRITECLUOUT, * )' PA3 .6, 10. 55;  ' 
WRITECLUOUT, * ) ’ PA . 3 , 10 . 55; ' 
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WRITECLUOUT, * D ' PA • 3 , 27 . 7; ’ 
WRITECLUOUT. *)'PA19.85.27.7; ' 
WRITECLUOUT. * ) ' PA19 . 85. 24 . 5; ' 
WRITECLUOUT. *D'PA24. 7, 24. 5;  ' 
WRITECLUOUT. * D ' PA24 . 7. 23 . 75;  ' 
WRITECLUOUT. *D'PA27. 15.23.75;  ' 
WRITECLUOUT. * ) ' PA27 . 15, 21 . 3 ;  ' 
WRITECLUOUT.  ’•D'PA27. 95, 21. 3;  ' 
WRITECLUOUT, * D ’ PA27 . 95, 8 . 0;  ' 
WRITECLUOUT, *  3 ' PA26 . 25, 8 . 0;  ' 
WRITECLUOUT, *3' PA26. 25. 7. 2; ' 
WRITECLUOUT, *  3 ' PA25 . 4 , 7 . 2 ;  ' 
WRITECLUOUT, *3'  PA25. 4, 6. 35;  ' 
WRITECLUOUT. *  3 ' PA24 . 6 , 6 . 35; ' 
WRITECLUOUT, *  3 ’ PA24 . 6. 5. 55; ' 
WRITECLUOUT, *  3 ' PA23 . 75, 5. 55; ' 
WRITECLUOUT, ‘ 3 ' PA23 . 75, 3 . 9;  ' 
WRITECLUOUT, *3'PA22. 9. 3. 9; " 
WRITECLUOUT, *  3 ' PA22 . 9,  3 . 15;  ' 
WRITECLUOUT, *  3 ' PA22 . 15, 3 . 15; ' 
WRITECLUOUT. *  3 ' PA22 . 15. 2 . 35;  ’ 
WRITECLUOUT, *3'PA21. 35, 2. 35;  ' 
W’TTECLUOUT, *3'PA21.35. .7;  ' 
WRITECLUOUT, *3' PA20. 55, .7;  ' 
WRITECLUOUT, * 3 ' PA20 . 55, - . 1 ; ' 
WRITECLUOUT. *3' PA10.05,-. 1; ' 
WRITECLUOUT. ‘3'PA10.05,  .7;  ' 
WRITECLUOUT, *3'PA8. 4, .7; ' 
WRITECLUOUT, *  3 ' PA8 . 4, 1 . 6; ' 
WRITECLUOUT, *3' PA6.8, 1.6;  " 
WRITECLUOUT, *3'PA6.8, .7; ' 
WRITECLUOUT,  ’‘3'PA6.05,  .7;PU;  ' 
RETURN 
END 


APPENDIX  III:  PIREP  Questionnaire. 


A  questionnaire  such  as  this  was  completed  by 
permanent-party  pilots  after  each  mission  flown  at  NIC 
during  two  2-week  training  periods  between  27  February  and 
12  May,  1988.  On  the  reverse  side  of  the  questionnaire  was 
a  contour  map  of  the  Ft  Irwin  area.  Data  from  those  areas 
annotated  by  the  pilots  were  consolidated  and  used  to  modify 
the  LTI. 


ON  THE  MAP  ON  THE  REVERSE.  REPORT  THE  FOLLOWING:  (NAME  OR 
AIRCRAFT  ID  NOT  REQDIRED  -  BE  AS  SPECIFIC  AS  POSSIBLE) 

1.  LOCATION. 

2.  TIME  (LOCAL). 

3.  ALTITUDE  ABOVE  GROUND  LEVEL  (AGL). 

4.  TURBULENCE  INTENSITY.  USE  STANDARD  REPORTING  CATEGORIES: 

NoneCN) ,  Light (L) , Moderate (M) , Severe( S) , extreme (X) 

PLEASE  INCLUDE  NEGATIVE  (N)  REPORTS. 


5.  SEE  EXAMPLE  ABOVE. 


APPENDIX  IV;  LTI  Inputs  and  Displays 


Tables  list  measured  surface  parameters  at  BYS  and  five 
sensor  locations  throughout  NTC,  and  interpolated  850  mb 
parameters.  Figures  are  hard-copy  reproductions  of  video 
displays  generated  by  the  LTI  program.  On  the  figures, 
Isopleths  of  threshold  LTI  values  for  turbulence  categories 
CCategory  I  aircraft)  are  overlaid  on  200  m  terrain  contours 
and  an  outline  of  Ft  Irwin.  Horizontal  ''tic"  marks  are 
drawn  at  2  km  Intervals. 
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Case 

Site 

Press 

Temp 

Dir 

Speed 

Cmb) 

C“F) 

Ckts) 

29  Feb  1200Z 

BYS 

899.9 

49 

267 

11 

APP  CmbD 

-0.5 

Tower  1 

921.5 

52 

285 

11 

850  Helaht  Cm)  1510 

Tower  2 

899.9 

49 

267 

11 

850  Temp  C“CD 

6 

Tower  3 

894.4 

49 

299 

4 

850  Wind  Dir. 

160 

Tower  4 

882.6 

48 

246 

9 

850  Wind  Spd  Ck.tsD 

15 

Tower  5 

937.7 

55 

100 

2 

01  Mar  OOOOZ 

BYS 

898.4 

56 

236 

12 

APP  CmbD 

-1.0 

Tower  1 

919.5 

61 

236 

19 

850  Height  Cm) 

1495 

Tower  2 

898.4 

56 

236 

12 

850  Temp  C“CD 

10 

Tower  3 

892.2 

58 

202 

16 

850  Wind  Dir. 

210 

Tower  4 

881.0 

54 

220 

13 

850  Wind  Spd  Ckta) 

20 

Tower  5 

934.6 

64 

227 

28 

01  Mar  1200Z 

BYS 

897.9 

44 

238 

6 

APP  CmbD 

-1.2 

Tower  1 

919.2 

48 

242 

12 

850  Height  Cm) 

1488 

Tower  2 

897.9 

44 

238 

6 

850  Temp  C-C) 

5 

Tower  3 

891.4 

47 

256 

9 

850  Wind  Dir. 

245 

Tower  4 

879.9 

44 

227 

10 

850  Wind  Spd  Chta) 

15 

Tower  5 

935.7 

50 

206 

8 

02  Mar  OOOOZ 

BYS 

893.9 

54 

250 

17 

APP  Cmb) 

-2.4 

Tower  1 

915.2 

59 

236 

19 

850  Height  Cm) 

1443 

Tower  2 

893.9 

54 

250 

17 

850  Temp  C^C) 

8 

Tower  3 

888.2 

54 

238 

10 

850  Wind  Dir. 

212 

Tower  4 

876.5 

52 

215 

13 

850  Wind  Spd  Ck.ts) 

16 

Tower  5 

931.0 

61 

227 

21 

05  Mar  1200Z 

BYS 

902.5 

53 

258 

6 

APP  Cmb) 

-0.4 

Tower  1 

923.9 

57 

279 

3 

850  Height  Cm) 

1527 

Tower  2 

902.5 

53 

258 

6 

850  Temp  CO 

12 

Tower  3 

896.7 

51 

345 

2 

850  Wind  Dir. 

30 

Tower  4 

885.1 

50 

43 

3 

850  Wind  Spd  Chts) 

3 

Tower  5 

939.8 

56 

150 

4 

06  Mar  OOOOZ 

BYS 

902.4 

70 

153 

4 

APP  Cmb) 

-1.8 

Tower  1 

923.3 

74 

205 

8 

850  Height  Cm) 

1533 

Tower  2 

902.4 

70 

153 

4 

850  Temp  C“C) 

13 

Tower  3 

896.5 

68 

171 

7 

850  Wind  Dir. 

340 

Tower  4 

885.2 

66 

173 

6 

850  Wind  Spd  ChtSj 

1  5 

Tower  5 

939.0 

75 

215 

4 

132 


06  Mar  1200Z 

APP  CabD 
850  Height  CmD 
850  Temp  C"C) 

850  Wind  Dir. 

850  Wind  Spd  Chts^ 

07  Mar  OOOOZ 

APP  CnbD 
850  Height  Cm) 

850  Temp  C"C) 

850  Wind  Dir. 

850  Wind  Spd  Chta) 

22  Apr  1200Z 

APP  CmbD 
850  Height  Cm) 

850  Temp  C^C) 

850  Wind  Dir. 

850  Wind  Spd  Chts) 

23  Apr  OOOOZ 

APP  Cmb) 

850  Height  Cm) 

850  Temp  C“C) 

850  Wind  Dir. 

850  Wind  Spd  Chta) 

25  Apr  1200Z 

APP  Cmb) 

850  Height  Cm) 

850  Temp  C^C) 

850  Wind  Dir. 

850  Wind  Spd  Chta) 

26  Apr  OOOOZ 

APP  Cmb) 

850  Height  Cm) 

850  Temp  C^C) 

850  Wind  Dir. 

850  Wind  Spd  Ck-ta) 


BYS 

900.7 

52 

261 

9 

-1.3 

Tower 

1 

921.8 

56 

254 

10 

1518 

Tower 

2 

900.7 

52 

261 

9 

12 

Tower 

3 

894.7 

54 

124 

8 

250 

Tower 

4 

883.2 

52 

199 

4 

5 

Tower 

5 

937.9 

57 

198 

7 

BYS 

896.9 

65 

236 

11 

-1.3 

Tower 

1 

917.8 

69 

234 

13 

1492 

Tower 

2 

896.9 

65 

236 

11 

14 

Tower 

3 

890.8 

66 

288 

6 

280 

Tower 

4 

879.8 

64 

241 

6 

10 

Tower 

5 

933.3 

71 

224 

18 

BYS 

890.6 

44 

248 

8 

0.6 

Tower 

1 

911.6 

47 

239 

17 

1433 

Tower 

2 

890.6 

44 

248 

8 

4 

Tower 

3 

884.6 

43 

130 

6 

280 

Tower 

4 

872.6 

42 

223 

7 

10 

Tower 

5 

927.7 

49 

215 

20 

BYS 

892.9 

56 

227 

17 

-0.2 

Tower 

1 

913.8 

60 

227 

27 

1461 

Tower 

2 

892.9 

56 

227 

17 

6 

Tower 

3 

886.7 

57 

214 

14 

245 

Tower 

4 

875.0 

54 

210 

18 

10 

Tower 

5 

929.5 

63 

219 

24 

BYS 

899.7 

50 

266 

9 

-0.4 

Tower 

1 

920.9 

53 

231 

16 

1509 

Tower 

2 

899.7 

50 

266 

9 

12 

Tower 

3 

893.5 

52 

158 

8 

280 

Tower 

4 

882.3 

48 

32 

2 

10 

Tower 

5 

937.0 

55 

191 

8 

BYS 

898.4 

73 

101 

5 

-1.5 

Tower  1 

919.3 

74 

57 

5 

1515 

Tower  2 

898.4 

73 

101 

5 

16 

Tower  3 

892.5 

71 

285 

4 

290 

Tower  4 

881.4 

70 

13 

6 

3 

Tower  5 

934.8 

78 

54 

2 

133 


26  Apr  1200Z 

BYS 

897.7 

64 

328 

2 

APP  CmbD 

-0.7 

Tower  1 

918.7 

64 

55 

3 

850  Height  CmD 

1489 

Tower  2 

897.7 

64 

328 

2 

850  Temp  C"CD 

15 

Tower  3 

892.0 

58 

67 

4 

850  Wind  Dir. 

180 

Tower  4 

880.7 

56 

40 

2 

850  Wind  Spd  Ck.tsD 

2 

Tower  5 

934.2 

62 

105 

3 

28  Apr  1200Z 

BYS 

892.3 

53 

275 

22 

APP  CmbD 

-0.7 

Tower  1 

914.2 

59 

248 

14 

850  Height  Cm!) 

1465 

Tower  2 

892.3 

53 

275 

22 

850  Temp 

10 

Tower  3 

887.5 

57 

293 

5 

850  Wind  Dir. 

220 

Tower  4 

875.9 

56 

251 

11 

850  Wind  Spd  ChtsD 

10 

Tower  5 

929.6 

62 

208 

18 

10  May  1200Z 

BYS 

903.5 

67 

20 

4 

APP  Cmb) 

0.1 

Tower  1 

924.5 

68 

236 

2 

850  Height  Cm) 

1549 

Tower  2 

903 . 5 

67 

20 

4 

850  Temp  C^CD 

18 

Tower  3 

897.8 

63 

270 

3 

850  Wind  Dir. 

40 

Tower  4 

886.5 

60 

10 

3 

850  Wind  Spd  ChtsD 

10 

Tower  5 

939.9 

69 

162 

4 

11  May  OOOOZ 

BYS 

903.4 

86 

103 

3 

APP  CnbD 

-1.8 

Tower  1 

923.9 

90 

209 

3 

850  Height  CaD 

1559 

Tower  2 

903.4 

86 

103 

3 

850  Temp 

22 

Tower  3 

897.3 

84 

354 

6 

850  Wind  Dir. 

25 

Tower  4 

886.7 

83 

20 

5 

850  Wind  Spd  Chta) 

10 

Tower  5 

938.9 

91 

348 

5 

11  May  1200Z 

BYS 

902.9 

74 

343 

9 

APP  Cmb) 

-0.2 

Tower  1 

923.7 

76 

61 

6 

850  Height  Cm) 

1555 

Tower  2 

902.9 

74 

343 

9 

850  Temp  C"C) 

20 

Tower  3 

897.4 

70 

44 

5 

850  Wind  Dir. 

350 

Tower  4 

886.2 

72 

320 

10 

850  Wind  Spd  Ckta) 

2 

Tower  5 

939.1 

72 

74 

4 

12  May  OOOOZ 

BYS 

901.0 

92 

94 

2 

APP  Cmb) 

-1.7 

Tower  1 

921.4 

93 

4 

3 

850  Height  Cm) 

1544 

Tower  2 

901.0 

92 

94 

2 

850  Temp  C“C) 

23 

Tower  3 

895.1 

88 

171 

4 

850  Wind  Dir. 

120 

Tower  4 

884.7 

88 

283 

4 

850  Wind  Spd  Chta]) 

2 

Tower  5 

936.1 

97 

325 

5 

134 


12  Hay  1200Z 

BYS 

900.3 

74 

74 

4 

APP  C«b) 

-0.3 

Tower  1 

920.9 

75 

294 

8 

850  Height  Cb^ 

1530 

Tower  2 

900.3 

74 

74 

4 

850  Temp  C“C) 

23 

Tower  3 

894.7 

74 

11 

3 

850  Wind  Dir. 

0 

Tower  4 

883.8 

73 

303 

6 

850  Wind  Spd  Chts^ 

0 

Tower  5 

936.2 

77 

182 

8 

136 


137 


138 


2  Mar  1988 
0000  GMT 


139 


5  Mar  1988 
1200  GMT 


6  Mar  1988 
0000  GMT 


143 


22  Apr  1988 
1200  GMT 


23  Apr  1988 
0000  GMT 


146 


Terrain  Contours 


147 


26  Apr  1988 


148 


149 


10  May  1988 
1200  GMT 


11  May  1988 
0000  GMT 


Terrain  Contours 


151 


152 


12  May  1988 
0000  GMT 


153 


APPENDIX  V;  Observed  vs  modelled  wind  speeds  and  LTI 

values . 


Values  are  given  for  the  five  meteorological  sensor 
locations  at  Ft  Irwin  for  case  days  and  times  listed. 

Table  V-1  shows  observed  winds  at  5  m  compared  to  10  m  winds 
output  by  the  WOCSS  wind  Interpolation  scheme.  In 
Table  V-2,  "observed"  LTI  values  are  those  calculated  using 
the  observed  winds  at  5  m,  while  "modelled"  values  Indicate 
LTI  values  generated  by  the  LTI  program. 
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Table  V-1 


Date/Time  CGMTD 

29 

FEB/1200 

11 

01 

MAR/0000 

19 

01 

MAR/ 1200 

12 

05 

MAR/1200 

3 

06 

MAR/0000 

8 

06 

MAR/12Q0 

10 

07 

MAR/0000 

13 

22 

APR/1200 

17 

23 

APR/0000 

27 

26 

APR/0000 

5 

26 

APR/1200 

3 

28 

APR/ 1200 

14 

10 

MAY/ 1200 

2 

11 

MAY/0000 

3 

11 

MAY/1200 

6 

12 

MAY/0000 

3 

12 

MAY/ 1200 

8 

Observed  vs  CModelled)  Wind  Speed  Ck.ts) 


TOWER 


1 

2 

3 

4 

5 

C13) 

11 

C15D 

4 

C16D 

9 

C16D 

2 

C14D 

C24D 

12 

C28D 

16 

C24D 

13 

C31D 

28 

C25D 

C16D 

6 

C17D 

9 

C14D 

10 

C19D 

8 

C16D 

C2D 

6 

C3D 

2 

C3D 

3 

COD 

4 

C2D 

C9:) 

4 

C9D 

7 

C9D 

6 

CIOD 

4 

C9D 

C8D 

9 

C9D 

8 

C6D 

4 

CIOD 

7 

C7D 

C13) 

11 

C13D 

6 

C13D 

6 

C14D 

18 

C12D 

C12D 

8 

C12D 

6 

CUD 

7 

C14D 

20 

CUD 

C36D 

17 

C36D 

14 

C38D 

18 

C37D 

24 

C37D 

C5D 

5 

C5D 

4 

C5D 

6 

C5D 

2 

C5D 

C2D 

2 

C2D 

4 

C2D 

2 

COD 

3 

C25 

C133 

22 

C15D 

5 

C12D 

11 

C19D 

18 

C13D 

CUD 

4 

C13D 

3 

C12D 

3 

C15D 

4 

CUD 

C3D 

3 

C3D 

6 

C3D 

5 

C2D 

5 

C4D 

C4D 

9 

C6D 

5 

C6D 

10 

C6D 

4 

C4D 

C2D 

2 

C2D 

4 

C3D 

4 

C2D 

5 

C3D 

CID 

4 

CID 

3 

CID 

6 

COD 

8 

CID 
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Table  V-2  Observed  vs  CModelled)  LTI 

TOWER 


12  3  4  5 


Date/Time  CGMT) 

29  FEB/1200  53 

C40) 

47 

C43D 

43 

C48D 

56 

C51) 

41 

C40) 

01 

MAR/0000 

69 

C56D 

55 

C62D 

63 

C56) 

68 

C73) 

75 

C56D 

01 

MAR/1200 

59 

C46D 

47 

C46D 

54 

C41) 

62 

C54D 

52 

C46D 

05 

MAR/ 1200 

42 

C35D 

38 

C47D 

38 

C47) 

47 

C12) 

40 

C355 

06 

MAR/0000 

70 

C823 

60 

C74) 

67 

CS2:> 

73 

C91) 

63 

C82:> 

06 

MAR/1200 

57 

C68D 

50 

C68D 

53 

C53:) 

57 

C83) 

52 

C61) 

07 

MAR/ 0000 

69 

C78:) 

61 

C72D 

60 

C72) 

67 

C83) 

71 

C72D 

22 

APR/1200 

63 

C52:) 

48 

C48) 

49 

C44:> 

58 

C60:> 

63 

C48) 

23 

APR/0000 

74 

CIOOD 

58 

C97D 

58 

C1025 

70 

C105) 

68 

CIOOD 

25 

APR/1200 

54 

C48D 

41 

C45D 

43 

C45D 

45 

C62D 

43 

C41D 

26 

APR/0000 

65 

C76D 

59 

C64D 

62 

C76:) 

71 

C89) 

59 

C76D 

26 

APR/ 1200 

52 

C29D 

45 

C43D 

51 

C43D 

56 

C14D 

49 

C29) 

28 

APR/ 1200 

61 

C59D 

62 

C63) 

49 

C55:) 

63 

C84) 

62 

C59D 

10 

MAY/ 1200 

41 

C36:) 

37 

C43D 

40 

C43:) 

47 

C53) 

40 

C40) 

11 

MAY/OOOO 

67 

C52) 

61 

C39D 

68 

C39) 

74 

C525 

66 

C52) 

11 

MAY/1200 

49 

C34D 

46 

C51D 

45 

C5i:) 

58 

C59) 

44 

C42:) 

12 

MAY/OOOO 

70 

C39D 

63 

C19) 

69 

C39) 

76 

C19) 

69 

C39:) 

12 

MAY/1200 

51 

C39D 

41 

C19D 

44 

C39:) 

55 

C19) 

49 

C39) 

